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We present a detailed report of the method, setup, analysis and results of a precision measurement 
of the positive muon lifetime. The experiment was conducted at the Paul Scherrer Institute using 
a time-structured, nearly 100%-polarized, surface muon beam and a segmented, fast-timing, plastic 
scintillator array. The measurement employed two target arrangements; a magnetized ferromagnetic 
target with a ~4 kG internal magnetic field and a crystal quartz target in a 130 G external magnetic 
field. Approximately 1.6 x 10^^ positrons were accumulated and together the data yield a muon 
lifetime of Tp(MuLan) = 2196 980.3(2.2) ps (1.0 ppm), thirty times more precise than previous 
generations of lifetime experiments. The lifetime measurement yields the most accurate value of the 
Fermi constant GF(MuLan) = 1.166 378 7(6) x 10"^ GeV"^ (0.5 ppm). It also enables new precision 
studies of weak interactions via lifetime measurements of muonic atoms. 



PACS numbers: 06.20.Jr, 12.15. Ji, 13.35.Bv, 14.60.Ef 

I. INTRODUCTION 

Weak interactions are germane to topics from big-bang 
nucleosynthesis and stellar hydrogen burning, to archaeo- 
logical dating and medical imaging, as well as fundamen- 
tal nuclear and particle physics. Unlike electromagnetic 
and strong interactions, all particles participate in weak 
interactions, and every process appears consistent with 
a universal strength for the weak force. This strength is 
governed by the Fermi constant [T]. 

The value of the Fermi constant Gf is best deter- 
mined by the measurement of the muon lifetime for 
two reasons; the muon lifetime is well suited to preci- 
sion time measurements and the pure-leptonic character 
of muon decay permits an unambiguous theoretical in- 
terpretation. Until recently, the uncertainty in Gp was 
17 parts-per-million (ppm); it originated from a 9 ppm 
contribution from the existing measurements of the muon 
lifetime and a 15 ppm contribution from the knowledge 
of the theoretical relation between the muon lifetime and 
the Fermi constant. 

The theoretical uncertainty arose from a limited knowl- 
edge of the 2-loop QED corrections to the muon decay 
rate. In 1999, van Ritbergen and Stuart (vRS) evaluated 
these corrections, thus establishing the relation between 



Gp and to sub-ppm level [2H1]- Their work was the 
motivation for our part-per-million measurement of the 
positive muon lifetime. 



A. Relation between Gf and the muon lifetime 

The V-A current-current Fermi interaction is a conve- 
nient effective theory for the description of the charged- 
current weak interaction and the definition of the Fermi 
constant. Although non-renormalizable, it is well suited 
to the low-energy limit of the electroweak interaction and 
neatly separates the parameterization of the weak inter- 
action from corrections due to electromagnetic interac- 
tions. In particular, Gp encapsulates all the weak inter- 
action effects in the low-energy effective theory [2^ . 

Feynman diagrams for ordinary decay at tree- and 1- 
loop level in Fermi theory are shown in Fig. [T] At the 
part-per-million level all appreciable corrections to the 
tree-level diagram arise from QED radiative corrections. 
Most importantly the hadronic effects are virtually non- 
existent, their leading terms arising as 2-loop QED cor- 
rections that contribute only 0.2 ppm to 

Despite such simplifications the extraction of Gf from 

to sub-ppm level requires an exhaustive calculation of 
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FIG. 1: Feynnian diagrams for ordinary muon decay in 
Fermi theory. The upper graph shows the tree-level di- 
agram for the ordinary decay amplitude while the lower 
graph shows the 1-loop QED correction to the ordinary 
decay amplitude. 



first- and second-order QED radiative corrections. In 
vRS, the authors expressed the relation between the 
Fermi constant and the muon lifetime as a perturbative 
expansion 

;r = ^(i+i:A,<->j (1) 

where is the muon mass and ^ Ag'*^ describes the 
phase space (Aq^^^), first-order QED (Ag^^^), second- 
order QED (Ag*^^^), etc., theoretical corrections. 

With an uncertainty arising from the theoretical cor- 
rections of 0.14 ppm [?-^ and an uncertainty arising 
from the muon mass of 0.08 ppm P|, the overwhelming 
uncertainty in the Fermi constant became the experimen- 
tal determination of the muon lifetime. 



B. Relation between Gf and the Standard Model 

Once Gi? is determined from t„ using Eqn.m it can be 
related to the fundamental parameters of any high energy 



The original derivation of vRS was conducted in the massless 
electron limit. The authors assessed a theoretical uncertainty of 
0.3 ppm in Gp arising from estimated mass terms at 2-loop order 
and neglected Feynman diagrams at 3-loop order. In 2008, Pak 
and Czarnecki [5] extended the work of vRS to include the 2-loop 
order, electron mass terms. Their work reduced the theoretical 
uncertainty in Gp to 0.14 ppm. 




FIG. 2: Standard model tree-level diagram of ordinary 
muon decay showing the W-boson mediating the elec- 
troweak interaction between the weak leptonic currents. 



theory, e.g. the electroweak interaction in the Standard 
Model. The relation between the Fermi constant and 
the Standard Model can be written as a perturbative 
expansion of the electroweak interaction 

where 1 /M^ represents the low-energy limit of the tree- 
level propagator corresponding to W^-boson exchange 
(See Fig. [2]) and ^ incorporates all higher-order elec- 
troweak interaction corrections [7]. Note that Eqn.[2]may 
be recast in a form with only observable quantities and 
the higher-order corrections 

This relation between the Fermi constant Gp, fine struc- 
ture constant a, and weak boson masses Mw and Mz, 
is an important prediction of the Standard Model. 

C. Summary of previous measurements 

In 1999, at the time of the vRS publications, the PDG 
world average of the muon lifetime was = 2.19703 ± 
0.00004 [S]. This average was largely determined by 
three measurements: Giovanetti et al. [5], Bardin et al. 
[TU] and Balandin et al. [TT], all conducted more than 
twenty five years ago. The experiments all observed the 
decay positrons from stopped muons in various targets. 

The experiments of Giovanetti et al. and Balandin et 
al. both used a continuous beam. To avoid any mis- 
identification of parent muons with detected positrons, 
the beams were restricted to low rates, thereby limiting 
the accumulated statistics of decay positrons. 

The Bardin et al. experiment used a pulsed beam 
with the associated time structure of the Saclay electron 
linac. In principle, a pulsed beam technique can over- 
come the aforementioned rate limitations by first accu- 
mulating multiple muon stops during the beam-on period 
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and then detecting multiple decay positrons during the 
beam-off period. In practice, the duty cycle at Saclay was 
not ideal and therefore limited the accumulated statistics 
of decay positrons. 

Following the vRS publications, two new measure- 
ments of the muon lifetime were initiated at the Paul 
Scherrer Institute. The FAST experiment fTT used a 
finely-segmented active target to identify decay positrons 
with parent muons and thereby operate at somewhat 
higher continuous beam rates. Herein we describe the 
MuLan experiment [13] . which made a part-per- million 
measurement of the positive muon lifetime using a cus- 
tomized time-structured muon beam. Letters describ- 
ing our 2004 commissioning measurement and 2006- 
2007 production measurements were published in Refs. 

[nmi]. 




FIG. 3: The positron energy distribution N{E) and 
energy-dependent positron asymmetry A(E) in ordinary 
muon decay fi~^ — e'^i'ei^fi- Em is the end-point energy 
of the ordinary decay spectrum. 



II. RELEVANT FEATURES OF DECAYS AND 
INTERACTIONS OF MUONS 

A. Muon decay. 

The known decay modes of positive muons are ordi- 
nary muon decay /i^ — >■ e^i^eP^, radiative muon decay 
^+ e+VeV^'y {BR = 1.4±0.4x 10^^), and the rare de- 
cay Ai+ e+e+e-i/eP^ {BR = 3.4 ± 0.4 x lO'^). In this 
experiment the neutrinos from muon decay were com- 
pletely undetectable, the photons from radiative decay 
were barely detectable, and the detected positrons were 
overwhelmingly from ordinary decay. Although the setup 
did not measure the energy of the positrons, it did iden- 
tify their direction. 

Due to parity violation in the weak interaction, the 
positrons emitted in ordinary decay are distributed asym- 
metrically relative to the muon spin axis. The distribu- 
tion is given by [T] 

N{e,E) cx 1 + A{E) cose , (4) 

where 9 is the angle between the muon spin direction 
and the positron momentum direction and A{E) is the 
energy-dependent positron asymmetry. As shown in 
Fig. |3j the asymmetry A{E) varies with positron energy 
E according to 

A{E) = {2E - Era)/{3Era - 2E) , (5) 

where Em = 1^,^/2 is the end-point of the positron en- 
ergy spectrum. For high energy positrons the asymmetry 
approaches A — +1 and positrons are preferentially emit- 
ted in the muon spin direction. For low energy positrons 
the asymmetry approaches A = —1/3 and positrons are 
preferentially emitted opposite the muon spin direction. 
When integrated over the theoretical energy distribution 
N{E) of emitted positrons, the asymmetry is A = -1-1/3. 



B. Muon slowing-down in condensed matter 

In the experiment a nearly 100%-longitudinally polar- 
ized surface muon beam was stopped in various target 
materials. The dynamics of a muon slowing-down in con- 
densed matter has important series for the initial proper- 
ties of the stopped muon ensemble in the target material. 

On entering the target material the incident muons are 
slowed by ionization and excitation. In the final stages 
of slowing-down, the muons can undergo a series of elec- 
tron pick-up and stripping reactions and exist as either 
neutral muonium atoms or charged muon "ions". Im- 
portantly, the formation of muonium can cause the rapid 
depolarization of stopping muons via the large hyperfine 
interaction between the /x/e magnetic moments. 

The initial depolarization and magnetic specie^ of 
stopped muons are consequently dependent on the details 
of the muonium formation and its hydrogen-like chemi- 
cal reactions in the stopping material. Empirically, large 
populations of diamagnetic /x"*" states are generally asso- 
ciated with muon stops in metals. By comparison muo- 
nium populations are less common, although their exis- 
tence is well established in such materials as ice, quartz, 
silicon, and germanium. Any surviving polarization of 
stopped muon or muonium populations is potentially 
worrisome for lifetime measurements. 



C. Muon spin rotation in condensed matter 

Following the slowing-down process the spin vectors of 
stopped muons will precess and relax in the local mag- 
netic fields of the target material. This phenomenon 



^ The initial depolarization refers to depolarization during the 
slowing-down process. The magnetic species refer to diamag- 
netic /i+ states and paramagnetic fj.^e~ states with compara- 
tively small and large magnetic moments, respectively. 
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is known as muon spin rotation (/iSR)0 This time- 
dependent muon-ensemble polarization yields, through 
the angular correlation between the muon spin and the 
positron direction, a time-dependent decay-positron an- 
gular distribution. In turn, this effect is seen — when de- 
tecting positrons in specific directions — as a geometry- 
dependent modulation of the exponential decay curve by 
the ^SR signals. 

In discussing muon spin rotation it is helpful to dis- 
tinguish the cases of /iSR in a transverse magnetic field; 
i.e., a magnetic field perpendicular to the polarization 
axis, and /iSR in a longitudinal magnetic field; i.e., a 
field parallel to the polarization axis. Of course, in gen- 
eral, both effects will contribute to the time distributions 
of the decay positrons. 



1. Transverse field (TF) fiSR 

A transverse field B causes the muon/muonium mag- 
netic moment jl to undergo Larmor precession in the 
plane perpendicular to the field axis. The precession 
frequency is a; = fiB/2mc where m represents the 
muon/muonium mass. It yields frequencies of 13.6 kHz 
per Gauss for free diamagnetic muons and 1.39 MHz per 
Gauss for free paramagnetic muonium|^ 

In matter, the ensemble-averaged polarization will gen- 
erally relax. One source of relaxation is the spin dephas- 
ing of the individual muons in the differing local i?-fields 
at their stopping locations in the target material. An- 
other source is the spin interactions between the stopped 
muons and the surrounding electrons and crystal lattice. 

For positrons detected at an angle Ob relative to the 
B-field axis, the effects of TF ^SR yield the time dis- 
tribution depicted in Fig. |4] A common, semi-empirical 
representation of the time dependence of the TF /iSR 
signal is 



N{t) cx (1 + A2{9b) e-^l^^ sin {ujt + 0)) e"*/^ 



(6) 



where ^2(^3) is the amplitude of the /xSR signal, uj and 
are the frequency and the phase of the precession, and Ti 
is the time constant of the relaxation. Note the geometry- 
dependent amplitude ^2(^5) is maximum for positrons 
emitted perpendicular to the i?-field axis and zero for 
positrons emitted parallel to the i?-field axis. Also — as 
shown in Fig.|4] — the TF /^SR signals of positrons emitted 
in opposite directions have opposite phases. 



Alternatively, /iSR can refer to muon spin rotation, relaxation or 
resonance. 

The muonium ground state is a two-state system consisting of the 
F = \ triplet state and the _F = singlet state. The precession 
frequency of triplet atoms is 1.39 MHz per Gauss. Singlet atoms 
are rapidly depolarized by hyperfine oscillations and therefore 
unobservable by /xSR techniques. 
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FIG. 4: Depiction of typical time spectra from TF /xSR 
(upper panel) and LR /zSR (lower panel). The paired- 
solid lines represent the time distributions of outgoing 
positrons in geometrically opposite detectors. The single- 
dashed line represents a positron distribution without 
^SR effects. The TF ^SR spectrum shows both spin 
precession and spin relaxation whereas the LF /iSR spec- 
trum shows a relaxation signal but no precession signal. 
Note, these figures greatly exaggerate the observed effects 
of /iSR in this experiment. 



2. Longitudinal field (LF) [iSR 

Because the torque f = jl x B on the muon spin van- 
ishes for a longitudinal _B-field, LF /iSR exhibits no pre- 
cession signal. However, in matter the polarization will 
relax via the spin interactions with the surrounding elec- 
trons and the crystal lattice. Because relaxation via spin 
dephasing in differing B-fields is absent, the time con- 
stants for LF-relaxation are generally longer than TF- 
relaxation. 

For positrons detected at an angle 9b relative to the 
_B-field axis, the effects of LF /iSR yield the time dis- 
tribution depicted in Fig. [4] A common, semi-empirical 
representation of the time dependence of the LF /xSR 
signal is 



N{t) oc {1±Ai{0b) e-*/^i) e-*/"- 



(7) 



where Ai (Ob) and Ti are the amplitude and the time con- 
stant of the longitudinal relaxation. Note the geometry- 
dependent amplitude Ai{6b) is maximum for positrons 
emitted along the i3-field axis and zero for positrons emit- 
ted perpendicular to the B-Held axis (i.e. the reverse of 
the TF-case). Again — as shown in Fig. |4— the LF fiSK 
signals of positrons emitted in opposite directions have 
opposite signs. 
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III. EXPERIMENTAL SETUP 

The experiment was conducted at Paul Scherrer In- 
stitute (PSI) in Villigen, Switzerland, using a nearly 
100% longitudinally polarized, surface beam. In- 
coming muons were stopped in a solid target and outgo- 
ing positrons were detected by a finely-segmented, large- 
acceptance, scintillator array instrumented with fast- 
sampling waveform digitizers. 

One important feature of the experimental setup was 
the time structure of the muon beam. This time struc- 
ture or "fill cycle" consisted of a 5 /is-long, beam-on ac- 
cumulation period followed by 22 /xs-long, beam-off mea- 
surement period. The time structure was important in 
avoiding the need to associate decay positron with parent 
muons, a requirement that limited earlier experiments 
using continuous muon beams. Rather, the time of the 
decay positrons are associated with the beam-off transi- 
tion. 

Another important feature was a large reduction of 
the muon ensemble polarization via the spin precession 
during the accumulation period. Since muons arrive ran- 
domly during accumulation, their spins precess through 
different angles under the influence of the local B-field, 
thus dephasing the individual spins and reducing the en- 
semble polarization. 

Additionally, the positron detector comprised a sym- 
metric array of target-centered, geometrically-opposite, 
detector pairs. As /xSR signals have equal magn itude s 
but opposite signs in opposite detectors (See Sec. II C I, 
any remaining effects from /xSR signals are largely can- 
celed by detector symmetry. 



A. Muon beam 

The experiment used the ttES beamline at the labora- 
tory's 590 MeV proton cyclotron. The 1.4-1.8 mA pro- 
ton beam provided the intense secondary muon beam 
required for this high statistics experiment. The pro- 
ton beam has a continuous macroscopic time structure 
consisting of roughly 1 ns duration proton pulses at a 
50.63 MHz repetition rate. 

The 7rE3 secondary beamline viewed the E-station pro- 
duction target in the primary proton beamline. The sta- 
tion comprised a 40 or 60 mm-long graphite target, which 
was rotated at roughly 1 Hz to prevent damage by beam 
heating. The rotation resulted in a slight modulation of 
the muon intensity. 

The surface muons originate from at-rest — )• /i^f^ 
decays in the outer layer of the production target. The 
decays yield a back-to-back neutrino and positive muon 
with momenta 29.8 MeV/c. As a consequence of parity 
violation in the weak interaction, the surface muons were 
nearly 100% longitudinally polarized, with their spins 
anti-aligned to their momenta. 

A schematic of the ttEB beamline elements from the E 
target station to the muon stopping target is shown in 



Fig. [5j The beamline includes a fixed-element section, 
which views the E target station at 90 degrees to the 
through-going proton beam. Surface muons produced by 
pion decays were collected by a magnetic lens system 
and momentum-selected through a 60-degree, vertically 
oriented, dipole. Additional quadrupoles and two hori- 
zontal and vertical slit systems provide momentum- and 
flux-limiting apertures. A second opposing dipole and a 
quadrupole doublet were tuned to redirect the beam par- 
allel to a raised experimental floor, which is roughly 6 m 
above the proton beam. 

The customized beamline elements that follow trans- 
port muons without deflection through the uncharged 
electrostatic kicker and next through an E x B veloc- 
ity separator, which removes unwanted positrons. The 
beam is then focused by a quadrupole triplet through an 
x-y slit system and then refocused by a second triplet 
at the stopping target. The three slit systems were ad- 
justed to optimize the combination of beam-on intensity, 
momentum acceptance, and beam-off extinction. 



1. Electrostatic kicker 

A time structure was imposed on the continuous 
beam using a custom-built, fast switching, electrostatic 
kicker [ini [T7]. With the kicker high voltage off, the 
muons were transported straight along the beamline axis. 
With the kicker high voltage on, the muons were dis- 
placed vertically by 46 mm and deflected downward by 
45 mrad at the kicker exit. They were then focused onto 
the bottom edge of the third slit system. 

The kicker system includes two pairs of 75 cm-long 
by 20 cm-wide aluminum plates that were aligned in se- 
ries with a spacing of 5 cm. The plates were positioned 
symmetrically above and below the beamline axis with a 
12 cm gap. Collimators with apertures of 12.5x12 cm^ 
were located both upstream and downstream to prevent 
particles reaching regions either above or below the plate 
gap. The plates were housed in a 60-cm diameter, 2.0-m 
long, cylindrical vacuum tank with high- voltage, vacuum 
feed-throughs. 

When energized, a potential difference of 25 kV was 
applied between the top plates and the bottom plates 
(+12.5 kV was applied to the two upper plates and 
— 12.5 kV was applied to the two lower plates with a 
virtual ground at the kicker mid-plane). The switching 
of plate voltages between ground potential and ±12.5 
kV (and vice versa) was accomplished using series cir- 
cuits ("stacks") of seventeen fast-transition, high-power, 
MOSFETs with each MOSFET switching 735 V. Pairs 
of MOSFET stacks were operated together in push- 
pull mode with one stack pushing the plate voltage to 
±12.5 kV and the other stack pulling the plate voltage 
to ground potential. 

The kicker was typically operated at a cycle frequency 
of about 37 kHz with a 5 /iS beam-on period (plates 
grounded) and a 22 /xs beam-off period (plates energized) . 



Experiment 



6 




FIG. 5: Schematic view of the 7rE3 beamhne including the E target station, opposing vertical bending magnets (D), 
quadrupole focusing magnets (denoted Q), slit systems (S), electrostatic kicker, E x B velocity separator, and the 
location of the experiment. See text for details. 



The transition time between ground and ±12.5 kV was 
67 ns [m [17] . Measurements performed during 2006 run- 
ning indicated the ±12.5 kV plate high voltage was stable 
to better than 300 mV. Improved measurements during 
2007 running indicated a stability of better than 150 mV. 
For further details see Sec. IVD 81 



grammable gate array (FPGA) that combined individual 
hits into x — y positions and time stamps of detected 
particles (during the accumulation period the data was 
pre-scaled). The FPGA data was read out through a 
Struck Innovative Systeme GmbH SIS3600 32-bit VME 
event latch 118] by the data acquisition system. 



2. Beam monitor 

A planar, high-rate, multi-wire proportional chamber 
with associated readout electronics was used for both ex- 
tensive beam studies prior to data taking and for peri- 
odic beam monitoring during the data taking. For the 
beam studies, the chamber was positioned at the stop- 
ping target location with the positron detector rolled to 
a downstream location. For the beam monitoring, the 
chamber was positioned immediately after the 150 /im- 
mylar beampipe window at the downstream end of the 
positron detector. 

The beam monitor was designed to measure the x-y 
muon profile during the high instantaneous rates of the 
accumulation period and the low instantaneous rates of 
the measurement period. The chamber was assembled 
from two perpendicular anode wire planes interleaved be- 
tween three high-voltage cathode planes. Each anode 
plane consisted of 96 tungsten anode wires of 20 /im 
diameter and 1.0 mm wire spacing and each cathode 
plane consisted of 12.5 /im aluminized Mylar foils held at 
3.0 kV. A "fast" chamber gas of 30% isobutane and 70% 
tetrafluoromethane was used to handle the high rates. 
The beam monitor was about 98% efficient for surface 
muons and about 8% efficient for beam positrons. 

Adjacent pairs of anode wires were connected together 
into 2 mm-spaced wire-pairs and read out by six, 16- 
channel, amplifier-discriminator boards. The resulting 
logic signals were recorded by a custom-built, field pro- 



3. Beam characteristics 

The beam time structure is depicted in Fig. [6| The 
upper panel shows the muon rate over the fill cycle and 
indicates the 5 /is beam-on accumulation period, 22.0 /iS 
beam-off measurement period, and the fast transitions 
between the beam-on/off states. The lower panel shows 
the corresponding positron rate over the fill cycle, which 
increases monotonically during the accumulation period 
and decreases exponentially during the measurement pe- 
riod. 

The muon beam was operated at a central momen- 
tum of 28.8 MeV/c and a momentum bite of ~2.6% (cr). 
During the beam-on accumulation period, the instanta- 
neous flux of surface muons at the stopping target was 
about 8 MHz in a spot size approximately 21x10 mm^. 
This rate yielded about 40 stopped muons per accumula- 
tion period and 15 detected positrons per measurement 
period. The surface muons arrive at the stopping tar- 
get with their polarization tilted at a 6°-angle below the 
beamline axis. This tilt arises from spin precession in the 
transverse i3-field of the velocity separator. 

During the beam-off measurement period the instanta- 
neous rate of surface muons was reduced by an "extinc- 
tion" factor e ~ 900 and the spot size was increased to 
about 43x32 mm^. 
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FIG. 6: Plot of the muon arrival times (upper panel) 
and decay positron times (lower panel) that are produced 
by the pulsed beam technique. The fill cycle comprises 
the Ta — 5 ^is beam-on accumulation period and the 
Tm = 22 fis beam-off measurement period. 



B. Target arrangement 

For the 2006/07 data taking the 7rE3 beamline was ex- 
tended through the positron detector with the stopping 
target mounted in the vacuum pipe. The in- vacuum tar- 
get reduced the number of upstream muon stops com- 
pared to our 2004 commissioning experiment 141 ■ These 
stops were worrisome as their /iSR signals are not can- 
celed by the opposite detector sums. Such distortions 
were identified in Ref. [Tij. 

The beamline extension necessitated a constriction of 
the vacuum pipe from 35.5-cm diameter to 20.2-cm diam- 
eter at a location 52 cm upstream of the stopping target. 
A hinge enabled the target to be swung into the beam 
path for production running and out of the beam path for 
beam monitoring. A collar enabled the stopping target 
to be rotated azimuthally about the beam axis. 

Two combinations of stopping target and transverse 
_B-fields were employed to maximize the spin dephasing 
of stopped muons. The first strategy used the large in- 
ternal B-field of a magnetized ferromagnetic alloy tar- 
get (Arnokrome-III) to rapidly dephase its population of 
diamagnetic muons. The second strategy used a moder- 
ate external S-field with a non-magnetic quartz crystal 
(Si02) to rapidly dephase its population of paramagnetic 
muonium atoms. 

To further reduce the dangers of upstream stops, the 
beampipe wall was also lined with Arnokrome-III ferro- 
magnetic foil. The 0.3 mm-thick lining was arranged with 
its magnetization perpendicular to the beamline axis and 
was extended 67 cm upstream of the stopping target. 



The lining dephases the ensemble-averaged polarization 
of upstream stops. 



1. Arnokrome-III target setup 

The Arnokrome-III (AK-3) stopping target was used 
in the 2006 production run. AK-3 is a ferromagnetic al- 
loy consisting of about 30% Cr, 10% Co and 60% Fe and 
manufactured by Arnold Engineering Co. The disk- 
shaped AK-3 target had a diameter of 200 mm and thick- 
ness of 0.5 mm and was mounted on the hinged frame in 
the vacuum pipe at the geometrical center of the positron 
detector. A SRIM [20^ calculation of the AK-3 target 
stopping power for the incident 28.8 MeV/c fi~^ beam 
gave a mean range of 0.17 mm and a range straggling of 
0.02 mm. 

The AK-3 target was magnetized in the plane of the 
disk. The internal i?-field was approximately 0.4 T and 
oriented for production running with the magnetization 
axis to the beam left or the beam right {i.e. at 90 degrees 
to the muon polarization). In the ^0.4 T transverse B- 
field the diamagnetic /i+'s have a ~50 MHz precession 
frequency. The resulting dephasing during the accumu- 
lation period reduced the ensemble-averaged transverse 
polarization by roughly a factor of 700. 

For further details on yitSR effects in Arnokrome-III see 



Sec. VB and Ref. [H]. 



2. Quartz target setup 

The quartz (Si02) stopping target was used in the 2007 
production run. The target was a disk-shaped, artificially 
grown, single crystal with a 130-mm diameter, 2.0-mm 
thickness that was oriented with its Z-crystallographic 
axis perpendicular to the disk face. The crystal was 
purchased from Boston Piezo-Optics Inc. [22] and grown 
by so-called seed-free technology to minimize crystal im- 
perfections that possibly lead to reduced formation of 
muonium atoms. A SRIM calculation of the quartz 
target stopping power for the 28.8 MeV/c muon beam 
gave a mean range of 0.53 mm and a range straggling of 
0.03 mm. 

The 130-mm diameter quartz disk was the largest di- 
ameter single crystal that was commercially available. 
The remaining region between the 130-mm quartz disk 
outer diameter and the 200-mm beam pipe inner diame- 
ter was therefore covered by a magnetized AK-3 annulus 
to guarantee depolarization of off-axis muon stops. 

A 130-G transverse B field was used for the spin pre- 
cession of the muon stops in the quartz target. The field 
was produced by a Halbach arrangement [23 of twenty 
four, 2.5x2.5x2.5 cm'^, neodymium magnets mounted in 
an aluminum ring of inner radius 14 cm, outer radius 
19 cm and thickness 7.5 mm. 

Note the magnet was mounted to the detector assem- 
bly rather than the beam pipe. The mounting scheme 
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FIG. 7: Cutaway drawing of the experimental setup 
showing the vacuum beampipe (1), stopping target (2), 
Halbach magnet (3), scintillator array (4) and the beam 
monitor (5) at the downstream window of the beam line. 
Note that the Halbach magnet was not installed for the 
AK-3 data taking and the stopping target was rotated 
out of the beam path for the beam monitoring. 



allowed for both its rotation about the beam axis and 
its inclination relative to the vertical axis. The rotation 
degree-of-freedom enabled the direction of the transverse 
B field to be varied and the inclination degree-of-freedom 
enabled the introduction of a longitudinal B field at the 
target. This setup required a careful alignment to ensure 
the muon polarization axis and magnetic field axis were 
mutually perpendicular during production running. 

In the 130-G transverse magnetic field the precession 
frequency of muonium atoms was roughly 180 MHz and 
that of positive muons was roughly 1.8 MHz. The result- 
ing dephasing during the accumulation period reduced 
the ensemble-averaged transverse polarization by roughly 
a factor of 1000 for the large paramagnetic muonium pop- 
ulation and 25 for the diamagnetic /i"*" population. 

For further details on fiSK effects in crystal quartz see 
Sec.lVBland Refs. 



C. Positron Detector 

The positron detector is depicted in Fig. [Tj The detec- 
tor was a fast-timing, finely-segmented, large-acceptance, 
plastic scintillator array that was geometrically centered 
on the stopping target. The detector enabled an inclusive 
measurement of positron time and angular distributions 
from the muon decays in the stopping target. The detec- 
tor segmentation was important in minimizing positron 
pileup and the detector geometry was important in han- 



dling /iSR distortions. 

The detector was constructed of 170 triangle-shaped 
scintillator pairs arranged in a truncated icosahedral 
(soccer ball) geometry. Each pair comprised an inner 
plastic scintillator tile and outer plastic scintillator tile. 
The pairs were grouped into ten pentagonal enclosures 
(pent-houses) containing five tile-pairs and twenty hexag- 
onal enclosures (hex-houses) containing six tile-pairs; the 
ten pent-houses and twenty hex-houses together forming 
the detector geometry. The 7rE3 beamline traversed the 
positron detector through one upstream pentagonal hole 
and one downstream pentagonal hole in the detector ar- 
ray. 

For later reference a coordinate system with the z-axis 
parallel to the beam axis and the x(y)-axes in the hor- 
izontal (vertical) directions perpendicular to the beam 
axis is introduced. The location of detectors are spec- 
ified by their polar angle 9 measured from the z-axis 
and their azimuthal angle (j) measured in the xy-plane 
(See Fig. [t]). The detector arrangement yielded pairs of 
geometrically-opposite tiles that viewed the target with 
coordinates {9, (j)) and (180° -9,(j) + 180°). 

The triangular tiles were 3-mm thick BC-404 scintilla- 
tor that produced light signals with 0.7 ns rise times and 
1.8 ns fall times |26j. The tiles in hex- houses were equilat- 
eral triangles with base lengths of 15.0 cm for inner tiles 
and 14.0 cm for outer tiles. The tiles in pent-houses were 
isosceles triangles with base lengths of 15.0 cm (13.8 cm) 
and heights of 10.0 cm (9.5 cm) for inner (outer) tiles. 
The scintillator array had an inner radius of approxi- 
mately 40.5 cm and subtended a solid angle of approxi- 
mately 70% of 47r steradians. 

Two edges of each tile were covered with reflective 
tape while the remaining edge was glued to an adia- 
batic lightguide followed by a 29-mm diameter photo- 
muhipher tube (PMT). The light guides and PMTs— not 
shown in Fig. [7] — emerged radially from the geometrical 
center of the detector array. Two types of photomulti- 
piers were used, with 296 tiles (eighteen hex-houses and 
eight pent-houses) instrumented by Photonis XP 2982 
PMTs [23 a-nd 44 tiles (two hex-houses and two pent- 
houses) instrumented by Electron Tube 9143 PMTs [28] . 
The same brands of 29-mm PMTs were used to read 
out geometrically-opposite detector pairs. A clip-line was 
used to shorten the durations of the pulses from the pho- 
tomulipliers to full-widths at 20% maximum of 10 ns. 

An important requirement was the accurate centering 
of the detector array on the stopping distribution. The 
detector was supported by a steel platform with an ar- 
rangement of vertical posts and horizontal rails enabling 
the fine adjustment of the detector horizontal and ver- 
tical coordinates. Using this arrangement, the detector 
array and stopping distribution were centered to better 
than 2 mm in each direction. 
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1. Detector operation 

On average each minimally-ionizing positron (MIP) 
yielded a photomultiplier signal of ^70 photo-electrons. 
The initial settings of the high voltages for each photo- 
multiplier were performed by adjusting the MIP peak to 
correspond to one half of the 1.0 V dynamic range of the 
readout electronics. During data taking the gain stabil- 
ity was continuously monitored via the amplitude of the 
MIP peak and showed typical drifts of ^2 mV per day. 
To compensate for long-timescale gain drifts, the PMT 
high voltages were adjusted on several occasions during 
each running periodP] 

The detector was designed to minimize vulnerability 
to gain changes. The vast majority of detected positrons 
were minimally ionizing particles with well-defined en- 
ergy loss in scintillator tiles. By constructing the detec- 
tor with two layers of scintillator tiles, and designing the 
detector for efficient collection of scintillation light, the 
resulting coincident pulses from through-going positrons 
were well separated from background noise. 

Fig. |8] shows the distribution of hits versus the angle 
relative to the beamline axis. The dip at 90 degrees 
is a consequence of the shadow of the target disk, target 
flange and Halbach magnet (in quartz data taking) . The 
increase in counts either side of 90 degrees is attributed 
to scattering in the same components. 



D. Laser system 

A laser system was used to monitor the time stability 
of the positron detector during the measurement period. 
The LN203 nitrogen laser [521 generated 337 nm wave- 
length, 600-ps duration, UV pulses that were split and 
then distributed via a fiberoptic cable to 24 scintillator 
tiles and one reference photomultiplier. The reference 
PMT was situated outside the experimental area in a lo- 
cation that was well-shielded from the beam-related ra- 
diation and the kicker-related electrical noise. The PMT 
pulse shapes corresponding to decay positrons and laser 
signals were similar. 

In normal operations the laser system was run asyn- 
chronously with the fill cycle at a 33 Hz repetition rate. 
Typically the system was operated for one hour each 8- 
hour shift. 



E. Waveform digitizers 

The PMT signals from scintillator tiles were read out 
via 340 channels of 450 MHz sampling-rate, 8-bit reso- 
lution, VME-based waveform digitizers (WFDs). Data 



xio 

S 8000- 
6000- 
4000- 

2000- AK-3 

Q 



' I I 
I I 



30 60 90 120 150 180 

(deg.) 



xlO° 

w 

S 4000- 
3000- 
2000- 
1000- 



■ i 



I i 



Quartz 



0, 



^ Slow gain changes over the data taking period do not interfere 
with the positron time distribution during the measurement pe- 
riod. 



30 60 90 120 150 180 

(deg.) 



FIG. 8: Hit distribution for the 170 tile pairs versus the 
angle 9 relative to the beamline axis for the entire AK- 
3/quartz datasets. The upper panel is the hit distribu- 
tion for the AK-3 setup without the Halbach magnet and 
the lower panel is the hit distribution for the quartz setup 
with the Halbach magnet. The vertical groupings with 
common ^-values are tile pairs with different values of 
angle 0. 



from the custom-built WFDs was used to determine the 
times and amplitudes of the tile pulses. The latter were 
also used in the detailed analysis of pulse pileup and gain 
changes. 

The WFD modules were single-slot VME-64 boards. 
Each module had four analog channels with 0-1 V ana- 
log input ranges, one enable input, and one clock in- 
put. The four channels of each module were used to 
digitize the signals from one inner-outer tile pair and its 
geometrically-opposite tile pair. Each enable input re- 
ceived a logic signal synchronized to the measurement 
period and each clock input received a 450 MHz sinu- 
soidal signal distributed by the clock system. 

The complete setup involved 85 WFD modules in six 
VME crates. The crates were located in racks on the de- 
tector platform in the experimental area. Approximately 
12 meter- long, 50 Q. cables connected the PMT outputs 
to the WFD inputs. 

Each analog input was continuously digitized by inde- 
pendent analog-to-digital converters (ADCs) that sample 
at the external clock frequency. Two field programmable 
gate arrays (FPGAs) controlled the data flow from the 
four sampling ADCs to their associated 512 kB, FIFO 
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memories. When the enable input was present and the 
analog input exceeded a threshold, a time stamp, fill cy- 
cle number, and contiguous block of 24 ADC samples 
were written to memory. The digitized data in the WFD 
memories were then read out via the VME backplane. 
The FIFO-design permitted simultaneous read and write 
access to digitizer memory. 

The threshold trigger that initiated storage of ADC 
samples was derived from a high-bandwidth voltage com- 
parator at each analog input. Note, while the ADCs 
digitized at the external clock frequency /, the FPGA 
control electronics operated at a divided frequency //4. 
Consequently, the recorded blocks of 24 ADC sample had 
ranges of 4-7 pre-samples preceding the "trigger" sample. 

Interesting waveforms occasionally exceeded 24 ADC 
samples, for example when two pulses were close in time 
or a single pulse with a large amplitude was over thresh- 
old for a long time. Therefore, if the WFD input was 
over threshold during the last ADC sample of the 24- 
sample block, the WFD was re-triggered and recorded a 
continuation block of 24 additional samples. 

The thresholds and offsets of each WFD channel were 
individually programmable via the VME interface. The 
WFD analog thresholds were set as low as tolerable for 
data taking rates. The WFD analog offsets were set 
to permit the recording of under-shoots following large 
pulses. 

Due to variations in clock signal cable lengths and 
WFD board propagation times, the relative timing of the 
enable signal and the (//4)-FPGA clock was not identi- 
cal on all modules. This not only results in a time offset 
between the different modules, but also a time offset that 
steps back and forth by 4 clock ticks depending on the 
relative timing of the enable signal and the (//4)-FPGA 
clock. Because the relative timing of the two signals 
changes for every measurement period, the displacements 
between WFD modules also change. Consequently, ac- 
cumulated time histograms were actually superpositions 
of physical time distributions with varying proportions 
of the 4 clock tick displacement. Note that because no 
ambiguity exists between the four inputs of a single dig- 
itizer, the definition of inner-outer tile coincidences and 
/iSR cancellation by opposite tile sums were unaffected. 
Most importantly, the determination of the lifetime was 
not affected by the displacements. 



F. Clock system 

An Agilent E4400 function generator [HU] was used for 
the external clock inputs of the waveform digitizer mod- 
ules. The function generator derives its clock signal from 
a temperature-controlled crystal oscillator with a quoted 
short-term and long-term stability better than 0.1 ppm. 

The clock was located on the detector platform in the 
experimental area. Its signal was amplified and split six 
ways to generate one copy for each WFD crate. The 
six copies were then amplified and passively split sixteen 



ways to produce one copy of the clock signal for each 
WFD module. 

A variation of the clock timebase that was correlated 
with the measurement period could cause distortions of 
the lifetime determination. Such effects could arise from 
beating between the clocks defining the digitizer sam- 
pling and kicker transitions. To eliminate this possibility 
the digitizer sampling frequencies and kicker transition 
frequencies were derived from completely independent, 
free running, clock systems. 

The collaboration was blinded to the exact frequency 
of the clock until the completion of the data analy- 
sis. During data taking and data analysis the frequency 
was known by the collaboration to be 451.0 ± 0.2 MHz. 
The unblinding involved first a relative unblinding of 
the two (2006 and 2007) datasets into commonly-blinded 
clock ticks and then an absolute unblinding of the entire 
dataset into real time units. 



G. Data acquisition 

The data acquisition ,31j provided the read out and 
data logging of the waveform digitizers and the beam 
monitor as well as run control, data monitoring, etc. The 
acquisition was developed using the midas acquisition 
package p2| and the ROOT analysis package [33 on a 
parallel, layered array of networked Intel Xeon 2.6 GHz 
processors running Fedora Core 3. 

To handle the high data rates and avoid any deadtime- 
related distortions, the readout was organized into re- 
peating cycles of dcadtime-free intervals called "data 
segments." Each data segment comprised a complete, 
deadtime-free record of every over-threshold digitized 
pulse during 5000 consecutive fill cycles. Between each 
data segment a short deadtime of roughly 2-4 ms was 
imposed to complete the transactions that packaged the 
data and synchronized the readout. The readout scheme 
relied on the ability to read out and write to the WFD 
FIFOs simultaneously. 

An important component of the acquisition system was 
a programmable pulser unit (FPU) that synchronized the 
data segments and measurement periods. The FPU was 
constructed from a 25 MHz Xilinx FPGA. It generated a 
continuous stream of fill cycles for the electrostatic kicker 
and a gated stream of 5000 fill cycles for the waveform 
digitizers. The data acquisition system initiated the be- 
ginning of each 5000 fill cycle sequence via the PPU and 
the PPU reported the completion of each 5000 fill cycle 
sequence to the data acquisition system. 

The data acquisition design involved a frontend layer 
for data readout, a backend layer for event building and 
data storage, and a slow control layer for control and 
diagnostics. The frontend layer consisted of seven dual- 
core processors for the parallel readout of the six VME 
crates containing the waveform digitizers and a single 
VME crate for the beam monitor readout. Each VME 
crate was read out using a Struck Innovative Systeme 
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GmbH SIS3100/1100 PCI-to-VME interface PJ. The 
frontend programs incorporated lossless compression of 
raw data using Huffman coding [351 as implemented 
in ZLIB [37]. 

The backend layer received data fragments asyn- 
chronously from the frontend layer via a 1 Gb Ether- 
net network. The fragments were assembled into com- 
plete events that comprised all the recorded pulses dur- 
ing the data segment. These events were then written 
to a RAID- 10 disk array and finally migrated to 400 GB 
LT03 tapes. 

The slow control layer was responsible for operation 
and monitoring of various instruments such as high 
voltage supplies, beamline magnet supplies, tempera- 
ture monitors, and i?-field monitors. Slow control data 
was read out periodically with one copy written to the 
main datastream and another copy written to an online 
database. 

Lastly, an online analysis layer was responsible for en- 
suring data quality. The online analyzer received data 
"as available" to avoid introducing any unnecessary dead- 
times into data readout. As well as recording basic quan- 
tities such as run start/stop times and event statistics, 
the online analysis layer determined quantities such as 
detector gains and fitted lifetimes. These run-by-run 
quantities were stored in an online database that pro- 
vided a comprehensive record of data taking and enabled 
a straightforward monitoring of data integrity. 



IV. DATA PREPARATION 

The procedure to obtain the positron time histograms 
from digitized scintillator pulses involves first determin- 
ing the various parameters of individual tile pulses and 
then constructing the coincidences between inner-outer 
tile hits. Finally, the accumulated histograms of coinci- 
dence times are corrected for systematic effects including 
pulse pileup and gain changes. 

The data obtained from the 2006 production run (R06) 
using the AK-3 target has 1.1 x 10^^ decay positrons and 
the data obtained from the 2007 production run (R07) 
using the quartz target has 5.4 x 10^^ decay positrons. 
The R06(R07) data were recorded in runs of duration 
2(10) minutes and size 2(10) GBytes, each containing 
about 4 X 10^(2 x 10^) positrons. In total they yielded 
about 130 terabytes of raw data. 

As listed in Table [ij the R06/R07 production data 
are divided into subsets with different orientations of the 
magnetic field relative to the beam axis and different po- 
sitions of the geometrical center of the detector relative 
to the surveyed intersection of the beam axis with the 
target disk. These subsets are important for consistency 
checks and systematic studies. 

Shorter systematics datasets were also collected. Ex- 
amples of systematics datasets omitted in the final statis- 
tics are (i) magnet orientations that introduced large 
longitudinal B fields at the stopping target, (ii) detec- 



tor alignments that introduced large geometrical asym- 
metries between opposite tile pairs, and (iii) shorter ac- 
cumulation periods that decreased the spin dephasing. 
Short measurements taken every shift with low rates of 
laser pulses illuminating a sub-set of scintillator tiles are 
included in the final statistics. 

The analysis was conducted using the 88-terafiop ABE 
cluster and mass storage system (MSS) at the National 
Center for Supercomputing Applications (NCSA) [35] . 
The processing of digitized waveforms into time his- 
tograms took roughly 60,000 CPU-hours for each dataset. 

TABLE I: Summary of R06/R07 production datasets in- 
dicating the subsets with different magnetic field orienta- 
tions and detector alignments. The magnetic field orien- 
tation is denoted L for beam- left and R for beam-right. 
The detector position represents the different coordinates 
of the geometrical center of the detector relative to the 
surveyed intersection of the beam axis with the target 
disk. 
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Prior to analysis, a number of checks of data integrity 
were performed for occasional hardware or software fail- 
ures. The checks ensured no loss or corruption of data 
between the digitizer electronics and the analysis code. 
About 1.2% of runs were rejected by integrity checks. 



A. Pulse fitting 

The time and amplitude of pulses are determined by 
least square fits to the individual digitized waveforms of 
each scintillator pulse. The procedure involves fitting a 
relatively high-resolution standard waveform (0.022 ns 
sampling-interval) to a relatively low-resolution individ- 
ual waveform (2.2 ns sampling-interval) and it assumes 
the pulse shape is independent of pulse amplitude. 



1. Pulse templates 

The high-resolution standard waveforms (pulse tem- 
plates) are constructed by combining a large number 
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of 450 MHz sampling-rate, single-pulse, digitized wave- 
forms. One template is derived for each scintillator tile 
using single-pulse waveforms with peak amplitudes of 
120-220 ADC counts. 

The template construction relies on the independence 
of the positron arrival times and the digitizer clock phase, 
which implies a uniform time distribution of pulse max- 
ima across the sub-sample time intervals between con- 
secutive clock ticks (c.t.). The waveforms of individ- 
ual pulses thus correspond to displaced sets of samples 
from the same underlying pulse shape. This makes it 
possible to reconstruct a high-resolution pulse template 
from many time-adjusted low-resolution waveforms via 
the computation of the aforementioned sub-sample time 
interval for each individual waveform. 

The procedure for computing the sub-sample time in- 
tervals of individual waveforms follows Ref. ^39J. It in- 
volves first calculating a sub-sample pseudo-time t' for 
each digitized waveform 



t' 



-tan 



Sm—1 



Sr, 



(8) 



from the three consecutive ADC samples Sm-i, Sm and 
Sm+i centered on the maximum sample Sm- Then, us- 
ing the measured distribution of pseudo-times t' and the 
uniform distribution of true times t, the pseudo-time is 
mapped to true time. 

Finally, using the true times the individual waveforms 
are time-aligned and superimposed in the construction 
of the pulse template for each scintillator tile. Fig. [9] 
(upper panel) shows a representative pulse template for 
an individual scintillator tile. 



2. Fitting algorithm 

The pulse fitting algorithm fits pulse templates to in- 
dividual waveforms so as to minimize 



samples 



n templates 



(9) 



where Si are the N samples of the digitized waveform 
and fi{tj) are the n templates of the scintillator pulses. 
The free parameters in Eqn. [9] comprise the time tj and 
amplitude Aj of each scintillator pulse and an assumed 
constant pedestal, P. In minimizing the quantity D, 
all samples are equally weighted, except for underflows 
or overflows which were excluded from the computation. 
The equal weighting of ADC samples is discussed further 
in Sec. IVD2l 

Given that the datascts contain ^10^^ pulses, the al- 
gorithm for fitting is optimized for speed and robust- 
ness. As the probability of two pulses on one wave- 
form is roughly 1%, the algorithm first tries to fit a sin- 
gle template to each waveform. Because the fit-function 
P + Aifi(ti) is linear in Ai and P, the minimum of D 



is found by a one-dimensional search on the non-linear 
parameter ti, with the other parameters calculated from 
their partial derivatives at the fit minimum. This one- 
dimensional search is performed by parabolic minimiza- 
tion using Brent's Method [40 . 

If the quality of this fit is not acceptable (as indi- 
cated by a search for further pulses in the residuals) 
the algorithm then tries to improve the fit by adding 
pulses, removing pulses and using more sophisticated, 
time-consuming, x^-minimization procedures. In fitting 
a waveform to n>l pulses, the fitter first tries an n- 
dimensional search on the pulse times ti using the Minuit 
minimization code |41j with the remaining parameters 
calculated from the aforementioned partial derivatives. 
If the exceeds a maximum value, the fitter attempts a 
2n+l dimensional search over all fit parameters. The fit- 
ter continues to add pulses or remove pulses until either 
(i) the residuals indicate no evidence for further pulses 
or (ii) ten attempts to fit n>l pulses are completed. 

During fitting, the algorithm will reject any pulses with 
amplitudes Aj less than 35 ADC counts, times tj outside 
the waveform boundaries, or times tj within three sam- 
ples of other pulses. When no pulses are found, the pulse 
fitter calculates an amplitude-weighted time of all ADC 
samples. 

Sample fits to digitized "islands" containing a single 
pulse and two overlapping pulses are shown in Fig.[9j The 
digitized waveforms were overwhelmingly single pulse is- 
lands with fractional contributions of roughly 1 x 10~^ 
from two pulse islands and roughly 7 x 10""* from zero 
pulse islands. 
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FIG. 9: Examples of the digitized waveforms and the 
corresponding fit results for a single pulse (upper panel) 
and two overlapping pulses (lower panel) . The waveforms 
are indicated by the histogram and the fit function by the 
curve. 
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B. Hit definition 

Two cuts are applied before constructing the coinci- 
dences; one defining an unambiguous software threshold 
for pulses and another defining an unambiguous software 
deadtime between hits. The construction of coincidence 
datasets with different threshold and deadtime settings 
is important for studying such effects as pulse pileup and 
gain variations. 

A typical pulse amplitude spectrum for an individual 
scintillator tile is plotted in Fig. [Toj It shows a broad 
peak due to decay positrons at ~120 ADC units and a 
rising background due to noise pulses below ~50 ADC 
units. As shown in Fig. 



10 these distributions are fit 



to the convolution of a Landau function and a Gaussian 
peak plus a linear background term. The fit-function 
minimum defines the normal amplitude threshold Athr, 
corresponding to the valley between the low-energy noise 
and the positron MIP peak. 

This threshold is re-calculated every 100(10) runs in 
R06(R07) data taking in order to account for the ob- 
served, long-timescale drifts in the detector gains over 
the data-taking period. Additional datasets, with higher 
thresholds of Athr + 20, Athr + 30 and Athr + 40, are pro- 
cessed for studying the effects of gain variations during 
the measurement period. 
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FIG. 10: Representative pulse amplitude distribution for 
an individual scintillator tile. The histogram represents 
the measured data and the smooth curve the fit function 
(see text for details). The fit-function minimum defines 
the normal amplitude threshold Athr for each scintillator 
tile. 



The pileup of positrons within the instrumental dead- 
time distorts the measured time distribution. Specif- 
ically, when two positrons are separated by >3 c.t. 
(6.7 ns), the corresponding pulses are cleanly resolved 
by the fitting algorithm. However, when two positrons 
are separated by <3 c.t. (6.7 ns), the algorithm cannot 
reliably discriminate the overlapping pulsesj^ 




To explicitly define the minimum resolving-time be- 
tween neighboring pulses, a software deadtime following 
above-threshold hits is introduced. It rejects any addi- 
tional hits within an artificial deadtime (ADT). A de- 
tailed study of pileup distortions with different ADTs of 
5-68 c.t. (11 — 151 ns) was made (see Sec. IV CI). 

A typical inner-outer pair time difference spectrum is 
plotted in Fig.jll] It shows a sharp coincidence peak from 
through-going positrons with a time resolution of 0.55 ns 
(1 cr). The time offsets between inner-outer pairs, origi- 
nating from differences in PMT voltages, cable lengths, 
and light guide configurations, are determined by a Gaus- 
sian fit to each inner-outer time difference histogram. A 
run-by-run table of timing offsets is used to account for 
changes of the offsets during the data-taking. 

The individual tile hits are then sorted into datasets 
of inner-outer coincidences, inner singles, and outer sin- 
gles. The inner-outer coincidences for each tile pair are 
identified using a coincidence window of ±1 ADT The 
coincidence time is defined as the time corresponding to 
the later hit of the inner-outer pair. 

In the unlikely case that two combinations of coinci- 
dence pairs are found within ±1 ADT — for example when 
one outer hit falls between two inner hits with mutual 
separations of 1-2 ADT — the earlier combination of hits 
defines the inner-outer coincidence (this reduces the sen- 
sitivity to coincidences originating from photomultiplier 
"ringing" ) . Such circumstances yield one coincidence hit 
and one single hit. 
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FIG. 11: Representative time difference distribution for 
a single inner-outer pair. The histogram is the measured 
distribution and the curve is the Gaussian fit. The ob- 
served timing jitter between inner-outer pairs is typically 
0.55 ns (1 cr). 

Finally, the sorted datasets of coincidence hits and sin- 
gles hits are used to fill per-run, per tile-pair, histograms 
of inner singles, outer singles and inner-outer coincidence 
times. 



Subsequently, an ADT of > 6 c.t. (13.3 ns) is used to avoid any 
re-triggering on the trailing-edges of the scintillator pulses. 



Using the same coincidence window and artificial deadtime sim- 
plifies the pileup correction procedure. It, however, does increase 
the contribution from accidental coincidences between unrelated 
singles hits in inner-outer tile pairs. 
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C. Histogram corrections 

Before fitting tlie time liistograms to extract tlie muon 
lifetime, several corrections are applied. The largest cor- 
rections account for time-dependent distortions arising 
from piieup of the decay positrons and variations of the 
detector gains. In addition, the statistical uncertainties 
on histogram bin contents are corrected for multiple tile- 
pair hits originating from positron scattering, positron 
annihilation and cosmic-rays. 



1. Pulse-pileup correction 

In filling the histograms, a well-defined, fixed-length, 
artificial deadtime is applied so that subsequent detector 
hits within the ADT of an earlier detector hit are not 
used. Consequently, the measured time distribution cor- 
responding to a parent time distribution N{t) cx e^*/'^^' 
will be 

N'{t) cx e"*/^^ - + ... (10) 




2 4 6 8 10 12 14 16 18 20 

time (|J,s) 



FIG. 12: Representative time distributions for the trigger 
hits before the piieup correction (curve 1) and corrections 
accounting for leading-order piieup (curve 2), extended 
piieup (curve 3), accidental coincidences (curve 4) and 
second-order piieup (curve 5). These histograms corre- 
spond to an artificial deadtime ADT = 6 c.t. (13.3 ns). 



where the second term represents the leading-order 
piieup losses having a magnitude r that depends on the 
beam rate and the ADT. 

The procedure for correcting piieup takes advantage 
of the time structure of the incident beam. The piieup 
losses are statistically recovered using so-called "shadow 
windows" , by replacing the lost hits in each measurement 
period with measured hits at equivalent times in neigh- 
boring measurement periods. 

To illustrate the method, consider the case of leading- 
order piieup, where a second hit within the artificial dead- 
time is lost. To correct this loss, if a hit is observed at 
a time ti in a fill j (denoted the "trigger" hit), a hit 
is searched for in an interval ti — tj -|- ADT in the fill 
j + 1 (denoted the "shadow" hit). Adding the resulting 
histogram of shadow hit times to the original histogram 
of trigger hit times thereby statistically restores the hits 
lost to leading-order piieup. 

The leading-order piieup correction for ADT = 6 c.t. 
(13.3 ns) is depicted in Fig. 12 (curve 2). The correction 
ranges from roughly 10"'^ at the start of the measurement 
period to roughly 10^^ at the end of the measurement 
period. 

The leading-order correction is insufhcient to entirely 
account for piieup effects. As the shadow windows have 
a width equal to the ADT, they contain either zero hits 
or one hit. Consequently — at the level of the higher- 
order piieup — the leading order pile-up is over-corrected 
and the higher-order pile-up is under-corrected by the 
shadow window (as > 1 hits in the shadow window are 
treated as 1 hit in the shadow window). 

To correct second-order piieup, if a shadow hit is found 
in the shadow window applied to fill j + 1 , another shadow 
hit is searched for in a second shadow window applied to 



fill J + 20 

The second-order piieup contribution for ADT — 6 c.t. 
is also depicted in Fig. 12 (curve 5). Its contribution 
is roughly 10~^ at the start of the measurement period 
and is completely negligible by the end of the measure- 
ment period. For a parent time distribution e"*/"^, the 
second-order piieup time distribution has a time depen- 
dence e^'^*/'^ and an amplitude that increases quadrati- 
cally with artificial deadtime. 

The piieup corrections described above assume the ab- 
sence of the a ~ 0.5 ns timing jitter between the inner 
and outer tile hits. This timing jitter causes additional 
piieup losses of inner-outer coincidences when one hit is 
lost but the other hit is not (recall the ADT is applied 
before defining the coincidence) . This "extended" piieup 
is measured by searching for single tile hits within the 
jitter of the trailing edge of each shadow window. 

The extended piieup correction is plotted in Fig. [12] 
(curve 3). Its contribution is roughly 10~^ at the start of 
the measurement period and is completely negligible by 
the end of the measurement period. For a parent time 
distribution e^*/'^, the extended piieup time distribution 
has the time dependence e~^*/'^ and an amplitude that 
is ADT- independent (at least for ADT < r). 

Random coincidences between unrelated singles hits 
can cause accidental coincidence hits. They can originate 
from muon decays, environmental backgrounds, and elec- 



* In practice, the shadow windows are applied using a five-fill ring 
bufli'er, with fill 1 being corrected by fills 2 and 4, fill 2 being 
corrected by fills 3 and 5, and finally fill 5 being corrected by fills 
1 and 3. 
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tronic noise. A significant number of inner tile hits with- 
out outer tile hits occur because inner tiles are slightly 
larger than outer tiles. Overall about 27% of inner hits 
are singles and about 3% of outer hits are singles. 

A shadow window technique is also used to correct for 
accidental inner-outer coincidences. An inner singles hit 
at time ti in the trigger fill initiates a search for an outer 
singles hit at time U — ti+ADT in the shadow fill. 

The measured contribution of accidental coincidences 
is depicted in Fig. 12 (curve 4). It ranges from roughly 
10~^ at the start of the measurement period and is negli- 
gible by the end of the measurement period. For a parent 
singles time distribution e~*/'^, the accidental coincidence 
time distribution has a time dependence e"^*/"^ and an 
amplitude that increases linearly with artificial deadtime. 

A Monte Carlo (MC) simulation is performed to ver- 
ify the procedure for correcting pileup. The MC gener- 
ates a sequence of beam cycles with tile hits from parent 
time distributions that reproduce the observed rates of 
coincidence hits and singles hits from muon decays and 
background sources. The simulation also incorporates 
the measured timing jitter between inner-outer pairs. 

The simulated data are analyzed using the same pro- 
cedures as the measured data, i.e., applying the artifi- 
cial deadtime, forming the inner-outer coincidences, and 
accumulating the histograms of trigger hits and shadow 
hits. The simulated trigger hit time distributions are cor- 
rected with the simulated shadow hit time distributions 
and are found to reproduce the parent hit time distri- 
bution. In particular, the "true" lifetime in the absence 
of pileup effects, accidental coincidences and timing jit- 
ter, and corrected lifetime after the accounting for pileup 
effects, accidental coincidences and timing jitter, are in 
good agreement. The lifetimes derived for ADT ranges 
of 5—68 c.t. (11—151 ns) are in agreement to a precision 
±0.1 ppm (see Fig. pi). 



2. Gain variation correction 

When applying deadtimes, constructing coincidences 
and filling histograms, only hits with amplitudes that 
exceed the threshold Athr are used. If the detector gain 
changes over the measurement period, then the time his- 
togram will be distorted, either by additional hits rising 
above the amplitude cut or by additional hits falling be- 
low the amplitude cut0 

The detector gain versus measurement time is obtained 
by accumulating a sequence of pulse amplitude spectra 
corresponding to consecutive, 220 ns-wide, time windows 
over the entire measurement period. Because the ampli- 
tude distributions for decay positrons and cosmic rays are 
different, and the relative proportion of decay positrons 
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^ Because the threshold Aj;^^ is applied to pedestal-subtracted am- 
plitudes the time histogram is not distorted by pedestal varia- 
tions. 



FIG. 13: Fitted muon lifetime versus artificial deadtime 
obtained from an analysis of a Monte-Carlo simulation 
of pileup effects, accidental coincidences and timing jit- 
ter. The simulated data are analyzed using the same 
software and methods as measured data. The fitted life- 
time is ADT- independent to better than 0.1 ppm and in 
statistical agreement with the "true" lifetime depicted at 
the location ADT = 0. The error bars on the ADT > 
points correspond to the statistical variations allowed by 
the ADT correction. The error bar on the ADT = point 
corresponds to the statistical variation allowed between 
the "true" lifetime that omits the pulse pileup, accidental 
coincidences and timing jitter, and the ADT > lifetimes 
that correct for pulse pileup, accidental coincidences and 
timing jitter. 



and cosmic rays is time dependent, the positron distribu- 
tions are obtained from measured distributions by sub- 
tracting the cosmic-ray distributions. The subtraction 
uses the measured cosmic-ray distribution in the last 
220 ns time window of the measurement period. 

The positron amplitude distributions are then fit to a 
convolution of a Landau energy loss distribution and a 
Gaussian instrumental resolution to derive the Landau 
most-probable-value (MPV) versus measurement time. 
In the absence of instabilities the MPV would be inde- 
pendent of the measurement time. 

Fig. [M]— which plots MPV versus measurement time 
for scintillator tiles sub-divided by photomultiplier 
type — indicates a gain variation of roughly 5x10^** over 
the measurement period. One feature, which appears in- 
dependent of PMT type, is an oscillation of the gain over 
the first microseconds of the measurement period. An- 
other feature, which appears different for PMT types, is 
a drift of the gain over the entire duration of the mea- 
surement period. 

A gain variation during the measurement period will 
distort the extracted lifetime. The change in the counts 
exceeding the amplitude cut, AA'^, and the change in the 
MPV of the amplitude distribution, AAmpv, sltb related 

by 

AN = 2 Nthr AAmpv [Athr/AMPv) (11) 
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FIG. 14: Fractional MPV change versus the hit time in 
the measurement period. The solid circles correspond to 
the amplitude variation of tiles with Photonis PMTs and 
the open triangles correspond to the amplitude variation 
of tiles with Electron Tube PMTs. 



where Nthr is the counts in the histogram bin correspond- 
ing to the amplitude cut, and AAmpv i^thr/^MPv) is 
the inferred change in the amplitude cut via a scaling 
of the measured change in the Landau MPV. The factor 
of 2 in Eqn. [iT] arises as gain changes of both inner and 
outer tiles can alter the above-threshold coincidences. 

Using the measured MPV versus measurement time 
and Eqn. [TT] the time histograms are then corrected for 
gain variations. In one approach, the measured AAmpv 
distribution is converted to a hit variation AA'^ and used 
to subtract (add) the extra (lost) hits from gain changes. 
In another approach, the MPV distribution is first fit to 
an empirical function describing the gain changes. This 
function is then converted from MPV versus measure- 
ment time to hits versus measurement time and used 
to make the gain correction (this approach is helpful in 
understanding the sensitivity of to oscillations and 
drifts). The gain corrections are applied separately to 
Electron Tube and Photonis PMTs. 

Gain instabilities during the measurement period are 
identified with the photomultipliers and the readout elec- 
tronics. 



(i) The gain oscillations in Fig. 14 originate from elec- 
tronics instabilities due to fan-out modules that dis- 
tribute the measurement period start signals to all 
waveform digitizer modules. The oscillations are 
reproducible in bench tests of waveform digitizers 
and are observed in all channels of all waveform 
digitizer modules in all VME crates. 



(ii) The gain drift in Fig. 14 is dependent on PMT 
type. To understand this effect, amplitude dis- 
tributions corresponding to different time intervals 
between consecutive detector hits are examined. 
The time dependencies of the Landau MPV ver- 
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FIG. 15: Time dependencies of the Landau MPV versus 
the time interval after a preceding pulse for represen- 
tative groups of Electron Tube PMTs (upper panel) and 
Photonis PMTs (lower panel). The black data points cor- 
respond to time intervals between two hits in the same 
tile pair and the same fill cycle. The red data points cor- 
respond to time intervals between two hits in the same 
tile pair but two neighboring fill cycles. The red points 
show no evidence of longer term effects from gain varia- 
tions between neighboring fills. 



sus the elapsed time after the preceding pulse are 
plotted in Fig. [15] for each PMT-type. The Pho- 
tonis PMTs show a relatively short-timescale ef- 
fect of duration ^25 ns and fractional gain change 
~1%. The Electron Tube PMTs show a relatively 
long-timescale effect of duration ~250 ns and frac- 
tional gain change ~1%. Because rates change with 
measurement time, these analog pileup effects con- 
tribute to the gain drift during the measurement 
period. 

The pileup effect in item (ii) is distinct from the "dig- 
ital" pileup associated with the artificial deadtime that 
was discussed earlier in Sec. IV C 1[ Effect (ii) is not 



corrected by the shadow window technique. 



3. Hit multiplicity correction 

The final correction accounts for multiple coincidence 
hits arising from single parent events. For example, a 
through-going cosmic ray or scattered decay positron 
may yield >1 coincidence hits. Consequently, the num- 
ber of coincidence hits in the time histograms exceeds 
the number of independent events populating the time 
histograms. 

To correctly account for statistical uncertainties, the 
separate time histograms of multiplicity-one coincidences 
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and multiplicity- two coincidences are accumulated p°] 
These histograms are then fit to obtain the relative con- 
tribution of multiplicity-one and -two coincidences to the 
decay curve and the time-independent background (the 
two contributions have different origins and different mul- 
tiplicities) . Denoting the multiplicity one and two con- 
tributions as Ni and N2 (decay term) and Ci and C2 
(background term), the ratio of independent events to 
coincident hits versus measurement time is 



R = 



(TVi + jV2)e~*/^>- + (Ci + C2) 
(A^i + 2iV2)e-*/^" + (Ci + 2C2) 



(12) 



The ratio R decreases from 0.97 at early times (i.e., a 
~3% contribution of multiplicity-two hits) to 0.96 at late 
times (i.e., a ^4% contribution of multplicity-two hits). 
To account for events with multiplicities >1 the Pois- 
son uncertainties on bin contents are inflated by 1/ \fR. 
Note that this correction only affects the uncertainty on 
the lifetime and the of the fit — not the value of the 
lifetime. 



V. LIFETIME ANALYSIS 

The extraction of the positive muon lifetime from fits 
to the positron time histograms is described next. Sup- 
plementary studies of transverse field and longitudinal 
field /iSR in AK-3 and quartz and various checks on data 
integrity are presented. Systematic corrections and sys- 
tematic uncertainties that arise from /iSR effects, pulse 
pileup, gain changes, and time-dependent backgrounds 
are described. 



A. Fitting procedures in AK-3 and quartz 

The muon lifetime is obtained by fitting the time his- 
tograms of inner-outer coincidences. In principle, each 
pileup-corrected, gain-corrected tile-pair time distribu- 
tion can be fit to the elementary function. 



N{t) = 7Ve~*/^^' -I- C 



(13) 



where is the muon lifetime 



is a normalization con- 
stant and C represents the flat background. In practice, 
distortions arising from residual /iSR effects must be con- 
sidered. Separate fitting strategies are employed for AK- 
3 and quartz because the specific characteristics of their 
/iSR signals are different. Most important is the absence 
(presence) in AK-3 (quartz) of discernible jjSK effects in 
the individual tile-pair time histograms. 



Multiplicity-one coincidences are defined as one inner-outer hit 
in an 8 c.t. (17.7 ns) interval and multiplicity-two coincidences 
are defined as two inner-outer hits in an 8 c.t. (17.7 ns) interval. 
The efi^ects of higher multiplicities are negligible. 
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FIG. 16: Three-parameter fit to the entire AK-3 dataset. 
The upper panel shows both the measured data and the 
fit function and the lower panel shows the normalized 
residuals between the data points and the fit function. 



In all cases the benchmark fits use histograms prepared 
with an artificial deadtime 6 c.t. (13.3 ns) and the afore- 
mentioned correction procedures for pulse-pileup, gain 
changes and hit-multiplicities. The nominal fit start and 
stop times are tstart = 0.22 /is and tstop = 21.12 /xs, a 
range which begins and ends at a safe interval from the 
beam (kicker) transitions. However, the results do not 
include the extrapolation to ADT = discussed in Sec. 
IV D II below. 



1. Fit to the AK-3 data 

A simple fit method is used for AK-3 production data. 
The summed tile-pair time histogram is fit to Eqn. |13[ the 
lifetime result thereby relying on sufficient cancellation of 
/iSR effects by spin dephasing, target choice and detector 
geometry. 

The benchmark fit to the entire AK-3 dataset (R06- 
A + R06-B) with Eqn. 



13 is plotted in Fig. 16 The 



extracted lifetime is = 2 196 980.1±2.5 ps. Both the fit 
X^/dof of 1224/1185 and the distribution of the residuals 
between the data points and the fit function indicate a 
reasonable fit and no evidence of any distortions such 
as uncanceled /iSR effects, mis-corrected pulse pileup or 
time-dependent backgrounds. 

The lifetime results for the benchmark fits to the sep- 
arate R06-A and R06-B datasets of r„ = 2 196 979.2 ± 
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3.4 ps and = 2 196 981.3±3.8 ps are in good agreement 
(the two datasets have reversed B-fields). The x^/dof of 
the fits and the distributions of the residuals are also 
acceptable. 



2. Fit to the quartz 



A more detailed approach to LF/TF /iSR effects is 
used to extract the lifetime from quartz data. First, 
geometry-dependent effective lifetimes are extracted for 
each tile-pair from fits using a modified version of Eqn. 13 



accounting for transverse-field precession. Then, is ex- 
tracted from the effective lifetimes via a fit function that 
accounts for geometry-dependent longitudinal-field relax- 
ation. Together these two steps account for all observed 
features of TF/LF /xSR in quartz. 

In step one the 170 time histograms of all tile-pairs are 
fit to 



N{t) = Ne-*^^''''[1 + fit)] + C 



(14) 



where Tos, N and C are the effective lifetime, normaliza- 
tion constant and background amplitude, respectively. 
The additional term f(t) accounts for the TF /iSR signal 
according to 



f(t) = Asm0B P2e-*/^= sin(a;i 



(15) 



where P2, T2, uj and are the initial polarization, relax- 
ation constant, angular frequency and phase of the TF 
/iSR signal. The quantity A is the asymmetry parameter 
of the e'''-angular distribution and the angle 9b is the 
tile coordinate relative to the B-field axis. The product 
Asm 9b determines the geometry-dependent amplitude 
of the TF /iSR signal in each tile pair (see Sec. II C I. 

A 5-parameter fit using Eqn. [14] is employed to extract 
the 170 tile-pair effective lifetimes. The parameters Tcs, 
iV, P2, (t> smd C are varied. The parameters T2 and w are 
fixed to the weighted averages of their best-fit values from 
an initial 7-parameter fit to a sub-set of tile-pairs having 
large TF signals. The asymmetry parameter is always 
set aX A = +1/3. Any difference between this theoretical 
energy-integrated asymmetry and the detector efficiency- 
weighted asymmetry is subsumed into P2- 

In step two the distribution of effective lifetimes versus 
tile-pair coordinates {9b, (I^b) about the B-field axis is fit 
to 



Teff(6lB, </.b) = r^(l + 5[9b, <i>B)) 



(16) 



where is the true muon lifetime and 



Bi4>b) is the 

tile-coordinate-dependent lifetime shift due to LF /zSR 
effects. The scheme for incorporating LF /iSR, which 



03 




20 40 60 80 100 120 140 160 180 

©B (deg-) 

FIG. 17: Effective muon lifetime Toff versus tile coordi- 
nate 9b relative to the i?-field axis. The data points are 
the 170 individual effective lifetime measurements and 
the solid curve is the effective lifetime fit from Eqn. [14] 
The small irregularities in the solid curve arise from the 
(/)B-dependence of the asymmetry parameter A. 



holds when Ti/Pi > r^p^ yields 

'Pi 



6{9B,<j)B) = , ^ 



A cost 



(17) 



where Pi is the initial longitudinal polarization and Ti is 
the longitudinal relaxation constant. The factor A cos 9b 
incorporates the asymmetry parameter A of the e+- 
angular distribution and the tile coordinate 9b relative 
to the P-field axis. It determines the amplitude of the 
LF /iSR signal in each tile pair (see Sec. II C I. 

a 2-parameter fit using Eqn. 
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IS per- 
and the initial 



To extract r^, 
formed varying both the muon lifetime 
polarization Pi. The relaxation constant Ti is fixed at 
28 ^s (see Sec. 



Pi 

"best-fit 



value for t 



V B 3 1 , a choice that infiuences the 



le polarization Pi but not r^. The 
asymmetry parameter A in Eqn. [17] depends on the co- 
ordinates {9B,<f>B) of the tile pair. This occurs because 
positron absorption and scattering processes by interven- 
ing materials yield a geometry-dependent angular and 
energy positron distribution in tile pairs. Therefore tile- 
dependent values of A{9b, 4>b) obtained from a GEANT 
simulation of the experimental setup are used. The asym- 
metries range from -1-0.30 to -1-0. 40 [^ 

Note the tile coordinates {9b,4'b) are calculated rela- 
tive to the center of the muon stopping distribution not 



In the quartz production data the longitudinal polarization is 
Pi < 0.01 and longitudinal relaxation constant is Ti = 28 it 8 ^s 
and therefore the condition Ti/Pi 2> is fulfilled. For further 
details see Sec. |V B| 

In practice the ditterence in between using the GEANT calcu- 
lated asymmetries A{9b,4'b) and theoretical energy-integrated 
A = 1/3 is 0.11 ppm. 



19 



(0 








■(5 


40 






£2. 


















30 


- 






20 








10 











, , , 1 


, , , , 



0.6 0.7 Q.i 



0.9 



1.1 1.2 



1.3 1.4 



FIG. 18: The distribution of x^/dof values from the 
170 fits to the individual time histograms for the R07- 
A dataset. The fit function (Eqn. [l4| ) incorporates the 
TF nSR signal. 



the detector array. The position offsets {Sx, 6y, Sz) of 
the stopping distribution from the detector center are 
obtained from the measured distribution of the outgoing 
positrons in the tile array. The offsets are determined to 
a precision of ±2 mm. 

Figs. [17] and [l8] show for dataset R07-A the fit to 
the angular distribution of the 170 effective lifetimes 
Tcs{0b,4>b) and the results for the x^/dof values from the 



170 tile-pair time fits. The distribution of the x^/dof val- 
ues from the time fits has an average 1.0, which demon- 
strates that TF /iSR effects are correctly handled. The 
fit of the Teff-distribution has a x^/dof of 201/168, which 
demonstrates the LF ^SR effects are resonably handled. 
The x^ of fits to decay time and Tcs distributions for 
other quartz datasets are of similar quality (see Table [lH 
Table [IT] summarizes the parameters derived from the 
application of this method to the quartz datasets. There 
is good agreement between the r^-values extracted from 
the six datasets. 

Because the six datasets involve different magnet ori- 
entations and different detector positions, the fitted val- 
ues of the initial longitudinal and transverse polarizations 
may differj^ The small values of Pi result from the per- 
pendicular arrangement of the muon spin and the i?-field 
while the small values of P2 results from the spin dephas- 
ing of the muon stops. 



B. Supplemental information on /iSR effects 

In AK-3 the 3-parameter fit is successful because /xSR 
signals are weak and cancel in sums of opposite tile pairs. 
In quartz the 5-parameter time distribution fits and 2- 
parameter t^b distribution fit accounts for all observed 
TF/LF /iSR distortions. We describe here supplemental 
information that further establishes that the /xSR effects 
in AK-3 and quartz were handled correctly. 



In particular, Pi is very sensitive to the exact angle between the 
Halbach magnet field and the muon spin axis. 



l-iSR studies in AK-, 



For decay times t > 100 ns, no evidence for TF /iSR is 
observed in AK-3. The absence of TF signals is expected 
as a consequence of the different axes and the varying 
frequencies of the spin rotation in the individual AK-3 
magnetic domains. 

A tiny TF fiSK distortion, compatible with the fast 
relaxation of the transverse polarization of the last muon 
arrivals during the accumulation period, is seen at times 
t < 100 ns. The signal is consistent with results from /zSR 
studies on AK-3 by Morenzoni et al. and Scheuer- 
mann et al. [42j . Both studies found evidence of spin 
relaxation with time constants of ^14 ns and ^50 ns, 
but no indications of longer timescales. 

As shown in Fig. [T9j a tiny variation of the fitted life- 
time versus the detector angle 9 relative to the beam axis 
is observed. This effect is attributed to the slow relax- 
ation of the longitudinal polarization Pi . 

The longitudinal relaxation axis is approximately par- 
allel to the beam axis while the AK-3 magnetization axis 
is oriented perpendicular to the beam axis. We believe 
this originates from the partial alignment of the AK-3 
magnetic domains along the AK-3 magnetization axis. 
Consequently, while Pi accrues only small contributions 
from stops in domains oriented nearly perpendicular to 
the /i-spin axis, it accrues much larger contributions from 
stops in domains oriented more parallel to the /i-spin axis. 
The result is an ensemble-averaged longitudinal polariza- 
tion that is roughly parallel to the beam axis. 

The LF /iSR effect on the fitted lifetime versus detector 
coordinate 9 is roughly 10 ppni. Supplemental data was 
acquired that enhanced Pi. The circular-cross section 
target, oriented at 90° to the beam axis, was replaced by 
an elliptical-cross section target, oriented at 45° to the 
beam axis. This increased the component of the AK-3 
magnetization along the /i-spin axis. The elliptical-target 
data shows a four-fold increase in lifetime variation with 
detector coordinate 6. This is qualitatively consistent 
with the expected increase in the polarization Pij^ It 
supports our argument that the ^-dependence of the fit- 
ted lifetime originates from a longitudinal relaxation in 
the AK-3 target, and therefore will cancel in sums of op- 
posite tile pairs. 



A quantitative prediction of the longitudinal polarization in- 
crease for the elliptical target requires a detailed knowledge of 
the domain orientations. However, the observed four-fold in- 
crease in the 9-variation of the fitted lifetime is consistent with 
a simple model of populations of fully aligned domains and ran- 
domly aligned domains that reproduces the ratio between the 
known AK-3 remnant and saturation fields ;19| . 
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TABLE II: Results of the fitting method for the six quartz datasets. The entries for the transverse polarization P2 are 
weighted averages obtained from tile-pairs having large TF /j,SR signals. The entries for the longitudinal polarization 
Pi are individual values obtained from effective lifetime fits. The parameter T2 — 4.7 /is is fixed from fits to time 
histograms with large TF /iSR signals and the parameter Ti = 28 fj,s is fixed from fits to supplemental data with 
amplified LF /iSR signals. The offsets {Sx, Sy, Sz) are fixed from the geometrical distribution of hits over the positron 
detector. 



parameter data set 

R07-A R07-B R07-C R07-D R07-E R07-F 

P2 X 1000 2.54 ± 0.02 2.45 ± 0.02 2.35 ± 0.05 2.51 ± 0.07 2.28 ± 0.05 2.54 ± 0.06 
w (Rad/^fs) 11.4440 ±0.0061 11.5342 ± 0.0066 11.570 ±0.016 11.489 ±0.019 11.579 ±0.016 11.498 ±0.018 
Pi X 1000 1.3 ±0.2 1.6 ±0.2 1.3 ± 0.4 1.2 ±0.7 0.7 ±0.5 1.5 ±0.6 
Ti, (ps) 2196981.93 ± 5.86 2196980.94 ± 5.51 2196975.43 ± 15.11 2196967.24 ± 22.79 2196985.33 ± 15.66 2196991.11 ± 20.40 
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2. Transverse field fiSR studies in quartz 
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FIG. 19: Plots of the fitted lifetimes from the individual 
time histograms grouped by the detector angle 9 relative 
to the beam axis. The upper panel represents the AK-3 
production data with the magnetization axis at 90 de- 
grees to the beam axis. The lower panel represents the 
AK-3 supplemental data with the magnetization axis at 
45 degrees to the beam axis. The four-fold increase in 
the ^-dependence of the fitted lifetime in the lower plot 
refiects the increase in the longitudinal polarization for 
the 45 degree dataset. Note the vertical scales on the two 
plots are the same. 




Figure 20 (upper panel) shows the ratio 



Ruoit) 



Nuit)-NDit) 
Nu{t) + ND{t) 



(18) 



involving the normalized time distributions of the coinci- 
dence hits Nij/£,(t) in the upstream/downstream hemi- 
spheres of the positron detectorj^ This diagnostic 
is sensitive to time-dependent differences between up- 
stream/downstream hits and clearly shows a TF preces- 
sion signal with a frequency w ^ 11 /.ts"^ and a time 
constant T2 ^ b fxs. The value of uj agrees with the 
expected value for a diamagnetic fi'^ population in the 
130 G Halbach magnet field. 

The relative amplitude of the TF /iSR signal in the 
production data is only ^10^'^. To further study the 
TF /iSR signal, a supplemental dataset was acquired by 
shortening the accumulation period from Ta = 5.0 to 
0.15 /xs and thereby reducing roughly 20-fold the spin 
dephasing during muon accumulation. As depicted in 
Fig. 20 (lower panel), the ratio Ruoit) obtained with 
Ta = 0.15 us shows a large TF signal consistent with a 
twenty-fold enhancement of the production data signal. 
Fits using Eqn. 14 to the time histograms derived from 
the Ta = 0.15 /is dataset indicate a quartz diamagnetic 
fi~^ population of roughly 6-7%, in good agreement with 
published work [Ml[^. 



A subtraction of the time-independent background and a nor- 
malization of the summed histogram contents, is performed be- 
fore deriving the ratio Ruoit) from the Nu/£){t) spectra. 
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FIG. 20: Ratio Runit) of the normalized dif- 
ference between the time distributions in the up- 
stream/downstream hemispheres of the positron detec- 
tor. The upper panel corresponds to production data 
with an accumulation period Ta — 5.0 /iS and the lower 
panel corresponds to supplemental data with a shortened 
accumulation period Ta = 0.15 ^s (note the different 
scales of the vertical axes). The plots highlight the con- 
tribution of the TF /^SR signal in the quartz data. 
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FIG. 21: Ratio of hit time distributions from the beam- 
left hemisphere {Ri^{t)) and the beam- right hemisphere 
(i?R(t)) to the detector sum. The upper panel corre- 
sponds to production data with a longitudinal polariza- 
tion Pi ^ 0.01 and the lower panel corresponds to supple- 
mental data with a longitudinal polarization Pi ^ 0.25. 
The gradual slopes in Rh/Rit) are the effect of the LF 
/LtSR in the quartz data. 



3. Longitudinal field /iSR studies in quartz 



The two ratios 



RUt) = NL{t)/{NL{t) + Nnit)) 
Rnit) ^ NR{t)/{NLit) + NR{t)) 



(19) 



are plotted in Fig. 21 ^^{t) and Nii{t) correspond 
to the normalized time distributions of the coincidence 
hits in the beam-left and beam-right hemispheres of the 
positron detector, respectivelyp^ These rates are sen- 
sitive to time-dependent differences between beam-left 
and right hits and suggest a gradual decrease (increase) 
in i?L(t) {Rji{t)) consistent with a long-timescale relax- 
ation of the muon longitudinal polarization along the B- 
field axis. The time constant Ti of the relaxation is much 
longer than the lifetime. 



^® A (i) subtraction of the time-independent background and (ii) 
normalization of the summed histogram contents is performed 
before deriving the ratios RL/R{t) from the Nj^/jiit) spectra. 



The relative amplitude of the LF /xSR signal in the 
production data is only 10"**. Supplemental datasets 
were acquired with enhanced LF effects by re-orientating 
the Halbach magnet from a perpendicular alignment to 
the /It"'" spin axis (Pi < 0.01) to non-perpendicular align- 
ments (Pi ~ 0.05- 0.25). 

shows the ratios Pl/_r(0 fo'" the Pi 



Fig. 
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0.25 



dataset and clearly indicates a gradual decrease (in- 
crease) in RL{t) {RR{t)) consistent with a roughly 
twenty-fold enhancement of the LF /iSR signal. Fits of 
Eqn.[T6]to the time histogram from the Pi = 0.05 to 0.25 
datasets give a relaxation constant Ti = 28 ± 8 /xs. 



C. Consistency checks 

A number of checks are performed on the overall con- 
sistency of the lifetime results. They include comparisons 
of fit results for different time ranges, run groups, detec- 
tor positions and P-field orientations. 



Fig. 22 shows the fitted lifetime versus the fit start time 
from tstart = 0.1 to 5.0 /iS for the AK-3 and quartz data. 
The points indicate the lifetime results from individual 
fits and the solid curves their permissible deviations (Icr) 
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from the benchmark fit with a start time tstart = 0.2 fis. 
The deviations account for the correlations between the 
individual fits and the benchmark fit. No evidence is 
seen for either a start-time dependence or a stop-time 
dependence of the lifetime. This builds confidence in the 
handling of known time-dependent effects including /iSR 
signals, positron pileup and gain changes, as well as the 
absence of any unidentified time-dependent distortions. 
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FIG. 22: Muon lifetime versus fit start-time for the AK- 
3 data (top) and the quartz data (bottom). The data 
points and gray band indicate the lifetime results and 
corresponding uncertainties from individual fits. The 
solid envelopes indicate the permissible statistical devi- 
ations (fcr) from the benchmark fits with a start-time 
tstart = 0.2 IIS. The la deviations are calculated with 
appropriate accounting for dataset correlations between 
different fit ranges. 



Fig. 23 shows (t^ — T^)/cri, the normalized deviations of 



lifetime values for the individual runs from the bench- 
mark result Tfj_. The figure indicates the normalized de- 
viations form a Gaussian distribution of mean zero and 
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FIG. 23: The distribution (t^ — T^)/ai of the normalized 
deviations of the fitted lifetimes from the individual 
runs from the final result r^. The data points are the 
run-by-run results and the solid curve is a least squares 
fit to a Gaussian distribution. The left panel is the AK-3 
dataset and the right panel is the quartz dataset. 



standard deviation 1.0, i.e., consistent with purely sta- 
tistical variations. These results for individual run varia- 
tions, and similar results for run group variations, demon- 
strate the absence of effects on from long-timescale 
changes in gains, pedestals and thresholds. 

Recall in AK-3 (quartz) tiny variations of the fitted 
lifetime with the coordinate 9 {9b) are observed. These 
effects are attributed to the slow relaxation of the lon- 
gitudinal polarization and are expected to cancel out 
in opposite-pair sums. However, such /iSR cancellation 
grows increasingly imperfect for a source that is displaced 
from the center of the detector. 

The 0(0B)-distribution of the lifetime values obtained 
by fitting summed time histograms of geometrically- 
opposite tile-pairs in AK-3 (quartz) are plotted in Fig. 
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Neither show any evidence for a systematic 9{9b)- 
dependence of the fitted lifetimes (as could originate from 
imperfect cancellation of /iSR distortions from upstream 
stops). 
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FIG. 24: The upper panel shows the lifetime results from 
opposite tile-pair sums grouped by detector coordinate 
9 for the entire AK-3 dataset. The lower panel shows 
the lifetime results from opposite tile-pair sums grouped 
by detector coordinate 9b for the entire quartz dataset. 
The plots demonstrate the cancellation of /^SR distor- 
tions from target stops in opposite tile pairs. 

A number of datasets were accumulated with different 
magnetic field orientations and different positron detec- 
tor offsets ((Sx, Jy, (5z). The agreement between lifetime 
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FIG. 25: Fitted lifetime for quartz datasets with different 
magnetic field orientations and positron detector offsets. 
Entries R07-A and R07-B correspond to the quartz pro- 
duction data with the i?-field orientation to the beam- 
left and beam-right, respectively. Entries R07-C through 
R07-F correspond to quartz production data with dif- 
ferent detector positions and entries 9 — —9.0° through 
9 = 20.0° correspond to quartz systematics data with 
different magnet tilt angles. The different detector posi- 
tions and magnet angles change the magnitudes and axes 
of the LF/TF polarizations and their corresponding /zSR 
signals. 



results for different magnet and detector configurations 
is an important requirement for the overall handling of 
the /xSR distortions. 

We conducted the most exhaustive studies on the 
larger /xSR signals from the quartz target. The quartz 
studies spanned a wide range of precession axes (from 
horizontal to vertical) and longitudinal polarizations 
(from ^ 0.001 to ^ 0.25). The specific datasets were 

(i) Production quartz data accumulated with the Hal- 
bach magnet B field oriented to the beam-left and 
the beam-right. This change reverses the preces- 
sion direction. 

(ii) Supplemental quartz data accumulated with the 
Halbach magnet i?-field tilted at pitch angles up 
to 40° from the vertical. This changes both the 
axis of the TF precession and the magnitude of the 
LF relaxation. 

(iii) Both quartz production and supplemental data 
accumulated with different offsets {dx, 5y, 5z) be- 
tween the detector center and the stopping distri- 
bution (the changes in offset were <2 cm). This 
changes the magnitude and direction of the Hal- 
bach i3-field at the stopping location and therefore 
the frequency of the TF precession and the ampli- 
tude of the LF relaxation. 

Fig.[25]summarizcs the lifetime values derived from the 
left /right B-field orientations, vertical/tilted i?-field ori- 
entations and different target-detector offsets. The good 



agreement between lifetime values obtained with differ- 
ent TF/LF /iSR signals verifies that muon spin rotation 
in quartz was handled correctly. 

Finally, we remark that the AK-3 production datasets 
were also accumulated with the target magnetization ori- 
ented to the beam- left and the beam-right. 



D. Systematics studies 

This section addresses four classes of systematic uncer- 
tainties. One class of effect is associated with corrections 
for pulse pileup and gain changes in the histogram con- 
struction and another class is associated with /xSR effects 
in the histogram fits. A third class is related to possible 
time-dependent backgrounds such as muon stops in the 
upstream beampipe and a time- varying extinction of the 
muon beam. A final class is the overall stability of the 
time measurement. 



1. Artificial deadtime correction 

A correction is applied for the lost hits during the ar- 
tificial deadtime (see Sec. IVCl). The pileup correc- 
tion is performed for artificial deadtimes (ADTs) of 5 to 
68 c.t. (11 to 151 ns). By increasing the ADT, one de- 
creases the trigger-hit histogram contents and increases 
the shadow-hit histogram contents. Assuming the proce- 
dure correctly accounts for digital pileup, the lifetime de- 
rived by summing trigger hit and shadow-hit histograms 
is independent of the applied ADT. 

Fig. 26 shows versus ADT both before and af- 
ter the pileup correction for the R07 dataset. Before 



the correction, the lifetime increases by 10.5 ppm/ns 
(R06) and 10.6 ppm/ns (R07), while after the correc- 
tion, the lifetime increases by 0.008 ppm/ns (R06) and 
0.007 ppm/ns (R07), a roughly 1000-fold reduction in the 
ADT-dependence of the lifetime r^. 

A "perfect" correction of pileup would yield an ADT- 
independent lifetime. The tiny residual slope in r^t ver- 
sus ADT after correcting for pileup thereby implies the 
procedure is slightly under- estimating the artificial dead- 
time effects at the level of about lO^'^. One possible ex- 
planation is either the omission of a small, high-order, 
pileup term or the mis-accounting for a larger, known, 
pileup term. However, as discussed in Sec. |IV C H the 
procedure for accounting for digital pileup was verified 
by Monte Carlo simulation and indicated no evidence of 
any omitted terms or mis-accounted terms at the 10"^ 
level. 

Another possible explanation is a non-statistical varia- 
tion of the muon stops in the accumulation periods. Non- 
statistical variations lead to a pileup under-correction for 
"high rate" fills and a pileup over-correction for "low 
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FIG. 26: Plot of versus ADT before the pileup correc- 
tion (open circles) and after the pileup correction (filled 
circles). The statistical uncertainties associated with the 
pileup correction are substantially smaller than the sta- 
tistical uncertainty on the muon lifetime and not visible 
on this scale. The solid lines through data points are 
only to guide the eye. 



rate" fills that does not cancel out exact lyP^ Although 
sources of non-statistical variations are known — e.g. ion 
source fluctuations and production target rotation — their 
effects on muon stops are insufficient to explain the slope 
in versus ADT. 

Although no conclusive explanation for the ADT- 
dependence of the fitted lifetime is identified, the ob- 
served variation of versus ADT is clearly linear. 
Therefore the lifetimes derived from ADT = 6 c.t. 
(13.3 ns) fits are corrected by linear extrapolations to 
zero deadtime. The extrapolation uses the slope obtained 
from a straight-line fit to versus ADT from 5 to 68 c.t. 
(11-151 ns). The extrapolation yields a decrease in of 
0.11 ppm for the AK-3 data and 0.10 ppm for the quartz 
data. A conservative systematic uncertainty of ±0.2 ppm 
is assigned to the digital pileup correction procedure. 



2. Below-threshold pulses 

In the absence of pileup, below-threshold pulses can- 
not generate hits. Two pileup effects involving below- 
threshold pulses are important: 

(i) "Pile-on" involves the pileup of two below- 
threshold pulses. If the pulses overlap, the 
below-threshold pulses can produce a single above- 
threshold pulse. Pile-on is a rate dependent in- 
crease in hits that distorts the time distributions. 



Because unaccounted pile-on adds hits at high rates 
(early times), it decreases the fitted r^. 

(ii) "Pile-down" involves the pileup of one below- 
threshold pulse and one above-threshold pulse. If 
the pulses do not overlap, the below-threshold pulse 
can raise the pedestal and lower the amplitude of 
the fitted above-threshold pulse. Pile-down is a rate 
dependent decrease in hits that distorts the time 
distributions. Because unaccounted pile-down re- 
moves hits at high rates (early times), it increases 
the fitted r^. 

The effects of pile-on in promoting below-threshold 
pulses, and pile-down in demoting above-threshold 
pulses, are not ADT-dependent and are not corrected by 
the shadow window technique. The two effects cause op- 
posite shifts in r^. For a pulse fitting procedure that 
is linear in the amplitude A and the pedestal P and 
equally weights all ADC samples, the effects exactly can- 
cel. However, the exact cancellation is defeated if below- 
threshold pulses are first identified in the pulse fitting 
and then dismissed by the amplitude cut. 

A complete evaluation of pile-on and pile-down re- 
quires the measurement of the below-threshold amplitude 
distribution. To avoid biases from the hardware thresh- 
old of the waveform digitizers, the distribution was mea- 
sured by acquiring data using an extended, 64-sample 
data block; the region of samples from 32 to 64 being 
used to obtain an unbiased amplitude distribution. A 
relatively large number of pulses having amplitudes <25 
ADC counts is found that carry no time dependence. A 
relatively small number of pulses having amplitudes 25- 
35 ADC counts is found that do carry the muon lifetime. 

Using the measured amplitude distributions, the com- 
bined effects of pile-on and pile-down on are estimated 
by Monte Carlo simulation to be less than 0.2 ppm. Im- 
portantly, the combined effects of pile-on and pile-down 
are manifest in the time dependence of the amplitude 
spectrum during the measurement period, and therefore 
are compensated by the application of the gain correction 
to the time histograms (see Sec. IVC2|. 



The effects do not cancel as the number of under-corrected hits 
in high-rate fills exceeds the number of over-corrected hits in 
low-rate fills. 



3. Time pick-off stability 

Pileup can result in the loss of valid hits or the addition 
of extra hits. Pileup can additionally "pull" the time 
determination of tile pulses through effects including the 
absorption of a small, overlapping pulse into the trigger 
pulse fit or the interference of a large, preceding pulse 
with the trigger pulse fit. Such distortions of the time 
"pick-off" are rate-dependent and therefore can distort 

The determination of the time pick-off stability is con- 
ducted using the 24 tiles instrumented with the optical 
fibers from the pulsed-laser system. Approximately 400 
laser runs were accumulated in each production run. The 
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time difference 6t{t) between laser pulses in the scin- 
tillator tiles and the reference PMT is sensitive to any 
rate-dependent shifts in the time pick-ofF from the tile 
pulsesp^ 

To determine the time difference St(t) versus measure- 
ment time t, a series of time difference histograms are ac- 
cumulated corresponding to sequential 2.5 /xs-wide time 
windows on the laser pulse time during the measurement 
period. These histograms are then fit to a Gaussian dis- 
tribution with the centroids determining 5t{t). The fit- 
ted slopes of St{t) versus t for all tiles are consistent with 
zero. By combining results from all tiles, the average 
time pick-off was found to change by less than 0.25 ps 
over the measurement period. The limit implies a negli- 
gible distortion of . 



ticular the ~10 ppm variation of the fitted lifetime with 
the tile-pair angle 6 about the beam axis. The cancel- 
lation is not exact owing to differences in the detector 
acceptances (for example, variations in tile efRciencies, 
electron absorption in intervening materials, and the mis- 
centering of the stopping distribution). An acceptance 
variation between opposite tile pairs of several percent is 
estimated from the measured variations of the positron 
rates in the scintillator tiles. By combining the observed 
(^-dependence of the lifetime with the estimated varia- 
tion of the opposite tile-pair acceptances, a 0.1 ppm sys- 
tematic uncertainty from AK-3 /iSR effects is derived. 



6. Quartz ^SR 



4- Gain variations 

For the normal amplitude threshold Atf^., the lifetime 
shifts resulting from gain variations are 0.50 ± 0.25 ppm 
(R06) and 0.53 ± 0.25 ppm (R07). The uncertainties of 
0.25 ppm are conservative and based on variations be- 
tween alternative procedures for making the gain cor- 
rections (see Sec. IV C 2). The separate shifts for pairs 
instrumented with the Electron Tube PMTs (22 tile 
pairs) and Photonis PMTs (148 tile pairs) are -1-5.0 and 
—0.2 ppm for the R06 dataset and -1-5.3 and —0.1 ppm for 
the R07 dataset. The separate values of obtained from 
Electron Tube and Photonis PMTs — after the appropri- 
ate correction for the gain variation — are in reasonable 
agreement for the two datasets. 

To test the procedure for correcting for gain varia- 
tions, time spectra having amplitude thresholds Athr+'^O, 
+30 and -1-40 ADC units were obtained. Increasing the 
threshold increases the sensitivity to gain changes due 
to a higher count fraction in the threshold amplitude bin 
(see Eqn. [u]). The Athr + 20, +30 and +40 settings 
increase the lifetime distortions from gain variations by 
factors of roughly 10, 30 and 50, respectively. A small 
error in the gain correction procedure at the normal am- 
plitude threshold thus implies a large error in the gain 
correction procedure at the higher amplitude thresholds. 
Lifetimes obtained from a 10% fraction of the R06/R07 
datasets with the Athr + '^O, +30 and +40 settings — with 
50-fold differences in gain corrections — are in reasonable 
agreement. 



5. AK-3 fj,SR 

The application of the fitting method to the AK-3 data 
relies on geometrical cancellation of /xSR effects, in par- 



As described in Sec. |V A2[ the handling of /iSR effects 
in quartz requires specifying the beam-detector offset, B- 
field axis and several /iSR parameters. 

In fitting the time histograms of tile-pairs, the pre- 
cession frequency uj and relaxation constant T2 are fixed 
to their "best fit" values (see Fig. 27 1. The change in 



when oj and T2 are varied by one standard deviation 
from their "best-fit" values is adopted as the associated 
systematic uncertainties. The procedure yields uncer- 
tainties in of 0.003 ppm from fixing uj and 0.02 ppm 
from fixing T2. 

The fit to the angular distribution of the effective life- 
times requires the determination of the tile coordinates 
{0b,4'b) relative to the i?-field axis. A ±2° uncertainty 
is assigned to the knowledge of the i?-field axis based 
on the uncertainties in the mechanical alignment of the 
Halbach magnet and the _B-field non-uniformities at the 
target locationj^ The uncertainty in is obtained by 
correspondingly varying the i?-field axis in the fit proce- 
dure for the effective lifetimes (Sec. V A 2 ). The approach 



yields an associated uncertainty of 0.004 ppm in r^. 

Also required in the fitting procedure to the effective 
lifetimes arc the offsets {Sx,Sy,dz) between the centers 
of the positron detector and the stopping distribution. 
These offsets are determined from measured differences 
between positron rates in opposite tile-pairs. A ±2 mm 
uncertainty is assigned to the determination of the off- 
sets {Sx, dy, Sz), based on the estimated variations in the 
tile-pair detection efficiencies. The corresponding uncer- 
tainty in is obtained by varying each offset by ±2 mm 
in the fit procedure to the effective lifetimes. The ap- 
proach yields an associated uncertainty of 0.19 ppm in 
T^. The dominant contribution originates from the offset 
6x along the direction of the longitudinal relaxation. 

Finally, an uncertainty in is associated with the 
positron asymmetries A{9b,4>b) required in the effec- 
tive lifetime fit procedure. This uncertainty is based 



The scintillator tile rate falls during the measurement period 
while the reference PMT rate is constant during the measurement 
period. 



The estimate of a ±2° uncertainty in the _B-field axis is consis- 
tent with the observed longitudinal polarizations in the quartz 
production data. 
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TABLE III: Compilation of the contributions to the sys- 
tematic uncertainty in the Ufetime from the fitting 
procedure for the /xSR effects in the quartz data. 
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FIG. 27: Examples of precession frequency w (upper 
panel) and relaxation constant T2 (lower panel) versus 
tile coordinate Ob obtained from 7-parameter fits using 



Eqn. 14 The solid horizontal lines indicate the fixed 
"best fit" values of uj and T2 that are used in the fits of 
the tile-dependent effective lifetimes. 



on the difference between obtained with the tile- 
dependent, GEANT-calculated asymmetries and the tile- 
independent, energy-integrated theoretical asymmetry 
A= 1/3. The approach yields an associated uncertainty 
of 0.05 ppm in t^. 

Table IIIII summarizes the various contributions to the 
systematic uncertainty in associated with the fitting 
procedure for the quartz /iSR effects. The total uncer- 
tainty in from quartz /xSR effects is 0.2 ppm. 



7. Upstream muon stops 

A small fraction of muons stop on the inner walls of 
the beampipe upstream of the target disk. These stops 
include both incoming muons in the beam halo and back- 
scattered muons from the stopping target. Both tar- 



parameter 


5t^ (ppm) 


UJ 


0.003 


T2 


0.021 


{ObAb) 


0.004 


5x 


0.19 


5y 


0.001 


5z 


0.001 


A 


0.05 



Total uncertainty 0.20 



get disks were sufficiently thick to eliminate downstream 
stops. 

The concern with upstream stops is potential shifts of 
arising from /iSR. Unlike target stops, where /iSR sig- 
nals cancel in opposite tile-pair sums, for upstream stops, 
the /iSR signals do not cancel out. Therefore the pres- 
ence of /iSR signals from upstream stops could distort 

Recall that following the downstream beamline el- 
ements, the vacuum-pipe diameter was reduced from 
35.5 cm to 20.2 cm at a location 52 cm upstream of 
the stopping target. To minimize any /iSR effects, a 
67 cm long section of the pipe immediately upstream of 
the stopping target was lined with 0.3 mm-thick AK-3 
foil. To maximize dephasing, the AK-3 magnetization 
was oriented perpendicular to the beam axis. 

Several measurements were conducted to estimate the 
effects of upstream stops. Two regions of beampipe are 
particularly worrisome; (i) the location of the beampipe 
constriction and (ii) the section of the reduced-diameter 
beampipe. 

To estimate the stopping fraction in the annular face 
of the pipe constriction, a 2x3 cm^ cross section, 3-mm 
thick, plastic scintillator was mounted on a vacuum feed- 
through in the beampipe region that was immediately 
upstream of the pipe constriction. The setup determined 
the radial dependence of the muon rate along horizon- 
tal, vertical and 45°-angled axes perpendicular to the 
beam direction. The measurements yielded an estimate 
of (0.5 ± 0.1)% for the stopping fraction in the pipe con- 
striction. 

To estimate the stopping fraction in the 20.2 cm diam- 
eter beam pipe, a plastic scintillator telescope was used 
to detect the outgoing positrons from beampipe stops. 
The telescope consists of two pairs of separated scin- 
tillator tiles that viewed a ~10 cm diameter region of 
beampipe. The measurement yielded a stopping fraction 
of (0.10±0.05)% in the 20.2 cm diameter beampipe. 

The combined results of measurements (i) and (ii) indi- 
cate less than a 1.0% upstream stopping fraction. To de- 
termine the effects of upstream stops, a dedicated dataset 
of upstream stops was collected. The dataset comprised 
^2x10^" stops in magnetized AK-3 at twelve upstream 
locations of 5 to 80 cm from the detector center. The 
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lifetimes derived from this upstream stop data and the 
production data are in good agreement. When the up- 
stream stop data and the production data are added in 
proportions of 1:100 the Ufetime decreases by 0.1 ppm. 



8. Beam extinction stability 

A time-dependent beam extinction (e ^ 900) will cause 
a time-dependent positron background. A variation of 
extinction can arise from a systematic change in the 
kicker voltage during the measurement period. 

A limit on lifetime distortions from voltage instabili- 
ties is derived by combining: (i) a measurement of the 
kicker voltage V(t) over the measurement period, and 
(ii) a measurement of the beam extinction e(V) versus 
the kicker voltage, with (iii) a simulation of the effects of 
a time-dependent background on the lifetime determina- 
tion. 

The kicker voltage was measured with a high voltage 
probe connected to the kicker deflector plates. The probe 
was read out using a digital oscilloscope and the average 
voltage V{t) from many kicker cycles was fit to determine 
a limit on the voltage drift during the measurement pe- 
riod. Data were obtained on several occasions during 
both production runs and yield limits of AV < 300 mV 
(R06) and AV < 150 mV (R07)(3 

The variation in beam extinction with kicker voltage 
Ae/AV was measured during both production runs us- 
ing the beam monitor. Over the kicker voltage range 
23.5-25.0 kV the extinction increased linearly with in- 
creased voltage. Combining the limits on AV^ with the 
determinations of Ae/AV, the limits on changes in beam 
extinction from voltage instabilities of Ae < 0.20 (R06) 
and Ae < 0.06 (R07) are obtained. 

A simulation is used to estimate the effect of an extinc- 
tion variation Ae on the fitted lifetime r^. First, a time 
distribution is generated, consisting of an exponential de- 
cay curve with a time-dependent background due to ex- 
tinction variation. Then, the resulting distribution is fit 
to an exponential decay curve with a time-independent 
background. The difference At^ between the simulated 
lifetime and the fitted lifetime thereby determines the dis- 
tortion of the lifetime originating from an unaccounted 
variation of the beam extinction. The procedure yields 
upper limits of 0.20 ppm (R06) and 0.07 ppm (R07) on 
the lifetime distortions from the kicker voltage drifts. 

These limits are 20-40 times smaller than the limit 
quoted in our 2004 commissioning measurement. The 
improvements in limits were achieved through better sup- 
pression of voltage drifts, better measurements of voltage 
stabilities, and increased beam extinction. 



9. Master clock stability 

Finally, the master clock output was compared at sev- 
eral frequencies to a calibrated rubidium frequency stan- 
dard. The comparisons were made at the beginning and 
the end of each running period. The observed deviations 
were less than 0.03 ppm over the duration of the exper- 
iment. A systematic uncertainty of 0.03 ppm on is 
assigned to the time calibration of the master clock. 



VI. RESULTS 



A. Summary of uncertainties 



Table IV summarizes the statistical and systematic un- 
certainties that contribute to the muon lifetime measure- 
ments in the R06 and R07 running periods. The uncer- 
tainties associated with fiSK effects and kicker instabili- 
ties are considered to be independent contributions to the 
R06 and R07 systematics. The larger fiSR uncertainty 
in R07 data-taking reflects the larger /iSR distortions in 
the quartz data and the smaller kicker uncertainty in R07 
data-taking reflects the improved limits on voltage insta- 
bilities. 

By comparison, the uncertainties associated with pulse 
pileup, gain variations, upstream stops, time pick-off sta- 
bility and clock calibration are considered to be common 
contributions to the R06/R07 running periods. In each 
case the same methods are used to obtain estimates of 
systematic uncertainties. 

TABLE IV: Sources of systematic uncertainties on the 
muon lifetime measurements in the R06/R07 running pe- 
riods. The uncertainties listed in single-column format 
are common uncertainties and those listed in two-column 
format are uncorrelated uncertainties. The last two rows 
are the combined systematic uncertainties and the overall 
statistical uncertainties for the R06/R07 datasets. 



Uncertainty 


R06 R07 
(ppm) (ppm) 


Kicker stability 


0.20 0.07 


fiSK distortions 


0.10 0.20 


Pulse pileup 


0.20 


Gain variations 


0.25 


Upstream stops 


0.10 


Timing pick-off stability 


0.12 


Master clock calibration 


0.03 


Combined systematic uncertainty 
Statistical uncertainty 


0.42 0.42 
1.14 1.68 



The "worst-case" voltage drift is obtained by examining a variety 
of fits [e.g. linear, exponential) to the kicker voltage during the 
measurement period. The values quoted for AV are obtained 
from a linear fit. 



B. Results for muon lifetime, r^j 

The individual results for the muon lifetime from the 
fits to the R06 (AK-3) dataset and the R07 (quartz) 
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dataset, after the —0.11 ppm (R06) and —0.10 ppm (R07) 
adjustments for the ADT = extrapolation, are 

r^(R06) = 2 196 979.9 ± 2.5{stat) ± 0.9{syst) ps (20) 

and 

T^(R07) = 2 196 981.2 ± 3.7{stat) ± 0.9isyst) ps. (21) 

The combined result 

T^(MuLan) = 2 196 980.3±2.1(stoi)±0.7(s2/st) ps (22) 

is obtained from the weighted average of the individ- 
ual R06/R07 values with the appropriate accounting for 
the correlated uncertainties. The result corresponds to 
an overall uncertainty in the muon lifetime of 2.2 ps or 
1.0 ppm. 

Figure [28] and Table |V] summarizes the measurements 
of the positive muon lifetime over the past forty years. 
They include the MuLan R06 and R07 results, the Mu- 
Lan commissioning measurement (Chitwood) |14j . the 
FAST commissioning measurement (Barczyk) |12| . and 
the earlier work of Duclos et al. [13], Balandin et al. 
[nHH, Giovanetti et al. |9] and Bardin et al. [10]. The 
final result from the combined R06 and R07 datasets is 
a factor of 10(15) improvement over the MuLan(FAST) 
commissioning measurements and a factor of 30-150 im- 
provement over the earlier experiments. 

TABLE V: Compilation of post-1970 measurements of 
the positive muon lifetime. The first three entries are 
our final results from the combined R06-I-R07 datasets 
and the individual R06/R07 datasets. Barczyk et al. is 
the result of the FAST commissioning measurement and 
Chitwood et al. is the result of the MuLan commission- 
ing measurement. The remaining entries are previous 
generations of lifetime experiments. In column one the 
uncertainty represents the combined statistical and sys- 
tematic errors. 



Measured 


Reference 


Publication 


lifetime (^s) 




year 


2.196 9803 ± 0.000 0022 


R06+R07 




2.196 9799 ± 0.000 0027 


R06 




2.196 9812 ±0.000 0038 


R07 




2.197 083 ± 0.000 035 


Barczyk '12' 


2008 


2.197 013 ± 0.000 024 


Chitwood ll4j 


2007 


2.197 078 ± 0.000 073 


Bardin [10| 


1984 


2.196 95 ± 0.000 06 


Giovanetti [5] 


1984 


2.197 11 ± 0.000 08 


Balandin [TTl El] 


1974 


2.197 3 ±0.000 3 


Duclos [43] 


1973 



The final result for is in reasonable agreement with 
the earlier work of Duclos et al., Balandin et at, Gio- 
vanetti et al. and Bardin et al. While it is in reasonable 
agreement with the MuLan commissioning measurement 
(a difference of 1.3 a) it is in marginal disagreement with 
the FAST commissioning measurement (a difference of 



• Balandin - 1974 

» Giovanetti - 1984 

• Bardin -1984 

——* Chitwood - 2007 

• Barczyk - 2008 

i\/1uLan - R06 

• iMuLan - R07 

2.19690 2.19695 2.19700 2.19705 2.19710 2.19715 2.19720 

(l^S) 

FIG. 28: Plot of post-1970 measurements of the positive 
muon lifetime. The data points represent the individual 
measurements and the vertical band corresponds to the 
combined R06 and R07 result. The result of Duclos et 
al. is not shown in the figure. 

2.9 a). The weighted average of all results in Table [v] 
gives a lifetime — 2196 981.1 ± 2.2 ps with a chi- 
squared x^/dof ~ 2.7. The weighted average is domi- 
nated by our final R06+R07 result with the poor chi- 
squared originating from the disagreement between our 
final result and the FAST commissioning result. 

C. Result for Fermi constant, Gf 

We use the relation obtained by van Ritbergen and 
Stuart (vRS) [2}3] for the extraction of the Fermi con- 
stant Gp from a measurement of the muon lifetime 
T^. The vRS relation is derived using the V-A current- 
current Fermi interaction with QED corrections evalu- 
ated to 2-loop order. It yields 

/ 1927r3 1 

\^ij^iTA^(^^yTA^^^yTA^ ^^^^ 

where is the measured muon lifetime, m^j is the mea- 
sured muon mass, and Aq^"), Aq<-^^ and Ag(2) 

are theo- 
retical corrections with Ag^^^ accounting for the effects 
of the non-zero electron mass on the fi-decay phase space 
and Aq^^/^^ accounting for the contributions of the 1-/2- 
loop radiative corrections to the /i-decay amplitude. The 
effects of non-zero neutrino mass are completely negligi- 
ble. 

For the muon lifetime, the weighted average of the 
R06+R07 results, = 2 196 980.3 ± 2.2 ps, is used. The 
determination of Gp in units of GeV~^ from the measure- 
ment of in units of ps requires a unit conversion via 
Planck's constant, h. For Planck's constant, the value 
recommended by the 2010 CODATA committee of 
h = 6.582 119 28(15) x lO'^s GeV • s is used. 

For the muon mass, the recommended value of — 
105.658 371 5(35) MeV [5] is used. This value is derived 
from the combination of the measurements of the electron 
mass and the electron-to-muon mass ratio. 

Computing the theoretical corrections requires both 
the electron-to-muon mass ratio me/m^ and the fine 
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structure constant a{m^) at the momentum transfer 
q = of the ^-decay process. The recommended value 
of mjm^ = 4.836 33166(12) x lO'^ [5| is used. For 
the fine structure constant a{m^) the value a{m^) — 
1.0/135.902 660 087(44) is used. The value is obtained 
from Eqn. 4.13 of van Ritbergen and Stuart [2. us- 
ing the CODATA value of the fine structure constant 
a(0) = 1.0/137.035 999 074(44) [5] at zero-momentum 
transfer. 

The value for the phase space term Aq^-^^ = 
— 187.1 ppm is obtained from Eqn. 2.7 in Ref. [2] 
and the value for the 1-loop QED correction Ag^^^ = 
-4233.7 ppm is obtained from Eqn. 2.8 in Ref. 0. Note 
that Eqn. 2.8 of Ref. [2| incorporates the effects of the 
non-zero electron mass on the 1-loop QED correctionj^ 

The value for the 2-loop QED correction Aq'^^^ — 
-1-36.3 ppm is obtained by summing the individual contri- 
butions from purely-photonic loops (Eqn. 9 in Ref. [3]), 
electron loops (Eqn. 10 in Ref. [3 ), niuon loops (Eqn. 19 
in Ref. d]), tau loops (Eqn. 20 in Ref. [1]) and hadronic 
loops (Eqn. 16 in Ref. f4j). Additionally, a correction of 
—0.4 ppm, first evaluated by Pak and Czarnecki is in- 
cluded to account for the effects of the non-zero electron 
mass on the 2-loop QED correction. 



Using Eqn. 23 and the aforementioned values for the 



and the theoretical corrections Aq^^\ Aq'-^'> and 



Agf^^ we obtain 

Gi.(MuLan) = 1.166 378 7(6) x 10^^ GeV 



(0.5 ppm). 

(24) 



This result represents a thirty-fold improvement over the 
1999 PDG value obtained before the vRS theoretical 
work and the lifetime measurements pre-dating MuLan. 

The error in Gp of 0.5 ppm is dominated by the 
1.0 ppm uncertainty of the muon lifetime measurement 
with additional contributions of 0.08 ppm from the muon 
mass measurement and 0.14 ppm from the theoretical 
corrections in Eqn. [23] Note the theoretical uncertainty 
of 0.3 ppm originally quoted by van Ritbergen and Stu- 
art [5H1] arose from estimates of electron mass terms in 
the 2-loop QED corrections (0.17 ppm), hadronic contri- 
butions in the 2-loop QED corrections (0.02 ppm), and 
3- loop QED corrections (0.14 ppm). Subsequently, Pak 
and Czarnecki [5 extended the work of vRS and cal- 
culated the 2-loop order, electron mass terms. Conse- 
quently, the overwhelming theoretical uncertainty is now 
the 0.14 ppm contribution associated with 3-loop QED 
corrections. The uncertainties associated with h, nie/mf^ 
and Q;(m^) are all negligible. 

The determination of Gi^ from has additionally been 
discussed by Erler and Langacker in the PDG Reviews 
of Particle Physics [46]. Their articles first appeared in 
1992 with later updates as measurements improved and 



theory advanced. The authors maintained a formulation 
that differs slightly from vRS with their placement of the 
phase space correction F{p) as a multiplicative factor for 
the radiative corrections, e.g. 



1 



1927r3 



F{P) 



(25) 

where p = rul/rnj^ is the electron-muon mass ratio 
squared, a(m^) is the fine structure constant evalu- 
ated at the relevant momentum transfer, and Hi{p) 
and i?2(p) represent the 1- and 2-loop radiative correc- 
tions. Although the phase space correction F{p) is iden- 
tical to 1 -I- Aq'^^\ and the 1-loop radiative correction 
F{p)Hi{p)a{m^) /tt is identical to Ag^^^ (see Ref. [45]). 
the 2-loop radiative corrections F{p)H2{p)a'^{mf^)/n'^ 
and Ag*^^^ differ by F{p) — thus introducing a small in- 



consistency between Eqn. [23]of vRS and Eqn. 25 of Erler 
and Langacker. However, as F{p) = 0.99981295 is almost 
unity, the evaluation of Gp using Eqn. [23] and Eqn. [25] 
gives the same value to a precision 0.1 ppm]^ 

The above extraction of Gf from is based on the 
assumption of the V ~ A chiral structure of the charged- 
current interaction. If the assumption of the V — A char- 
acter is relinquished, the relation between and Gp is 
modified l47l. 



D. Discussion of results 

Below we compare the results from the AK-3 (R06) 
dataset and the quartz (R07) dataset and comment on 
the significance of the 1.0 ppm measurement of the muon 
lifetime and the 0.5 ppm determination of the Fermi con- 
stant. 



1. Comparison between R06/R07 datsets 

The two values of extracted from the R06/R07 
datasets are in good agreement. Because the two datasets 
involved different targets having greatly varying mag- 
nitudes, time dependencies and angular distributions of 
/iSR signals, the natural interpretation is that the spin 
rotation aspects have been properly handled. 

However, a second concern for a precision lifetime mea- 
surement is any difference between the positive muon 
lifetime in free space and the positive muon lifetime in 
matter. For example, a muon can reside in matter as 
either a bound muonium atom or an unbound positive 
muon. The relative populations of bound muonium and 



The 1-loop QED correction including electron mass effects was 
originally derived by Nir I45l . 



In addition, the expressions used by vRS and Erler and Lan- 
gacker to evaluate the running of a(m^) are slightly different. 
However, these differences are numerically unimportant for the 
Fermi constant at the part-per-million level. 
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unbound /x"'"'s can vary with material and ranges from 
small muonium fractions {e.g. AK-3) to large muonium 
fractions {e.g. quartz). A significant difference in life- 
times between muonium and /i"'"'s could jeopardize the 
interpretation of lifetime measurements using stopped 
muons. 

We stress: the two targets involve different populations 
of /x"'"'s (the dominant species in AK-3) and muonium 
atoms (the dominant species in quartz). The possibility 
of a lifetime difference between the free muon and the 
muonium atom was addressed by Czarnecki, LePage and 
Marciano ^S]- They found that the leading-order con- 
tributions from the two largest corrections — the /x+e" 
binding energy and e+e^ final-state interaction — cancel 
out and the resulting difference between their lifetimes 
to be only ~0.6 parts-per-billionj^ This difference is 
negligible compared to the part-per-million precision of 
the R06/R07 lifetime measurements. 

In addition, a lifetime difference between free muons 
and muonium atoms could arise from /i+e~ annihilation 
in muonium atoms. For example, the rate of ii^e~ — >■ 
v^Ve annihilation was considered by Czarnecki et al. |48| 
and earlier by Li et al. |51j . However, the SM prediction 
for the annihilation rate is only 10~^^ of the muon de- 
cay rate and therefore negligible compared to the ppm 
precision of the RG6/R07 lifetime measurements. 

2. Relevance to precision muon capture studies 

A precise determination of is important to current 
work on nuclear muon capture. The MuCap experiment 
|52j has recently measured the /x~p singlet capture rate 
As and the MuSun experiment [5 3) is presently mea- 
suring the /J,~d doublet capture rate A^. Both experi- 
ments derive the muon capture rates from the tiny dif- 
ference (Ao — A) between the positive muon decay rate 
(Aq = l/r^+) and the muonic atom disappearance rates 
(A = 1/t^i-z). Because the capture rates are very small, 
a precise determination of the muonic atom disappear- 
ance rates and the free muon decay rate is necessary. 

Prior to this work the muon decay rate Aq was known 
to 9 ppm or 5 s^^, an uncertainty that limited the cap- 
ture rate determinations from disappearance rate experi- 
ments. The new result for the muon lifetime — yielding a 
precision of 1 ppm or 0.5 in A^ — has eliminated this 
source of uncertainty. 

3, Relevance to fundamental constant determinations 

In 1999 the uncertainty in Gf was 17 ppm, a combina- 
tion of a 9 ppm experimental uncertainty and a 15 ppm 
theoretical uncertainty. The combined outcome of the 



vRS calculations [3] and the MuLan experiment is a 
thirty fold reduction in the total uncertainty in the Fermi 
constant — to a precision 0.5 ppm. 

The electroweak sector of the standard model involves 
three parameters, the two gauge coupling constants g, gi 
and the Higgs energy density v. Their values are fixed 
by measurements of the fine structure constant, a, Fermi 
coupling constant, Gi?, and Z boson mass, Mz- The 
thirty- fold improvement in the determination of Gf, to- 
gether with recent improvements in the determinations 
of a and Mz, allows for improved tests of the Standard 
Model. 



VII. CONCLUSION 

In summary, we report a measurement of the posi- 
tive muon lifetime to 1.0 ppm and a determination 
of the Fermi constant Gp to 0.5 ppm. The experiment 
used a novel, time-structured muon beam with a seg- 
mented, fast-timing, positron detector instrumented with 
450 MHz waveform digitizer readout. The results for 
and Gp are substantial improvements over earlier deter- 
minations from previous generations of lifetime experi- 
ments. Improved knowledge of is important to other 
experimental programs and improved knowledge oi Gp 
is important for numerous electroweak processes. 

In the determination of the Fermi constant, we are 
fortunate that the uncertainties on many quantities-e.g. 
the muon mass and Planck's constant — are significantly 
below the ppm-level. To improve the determination of 
Gp by another order of magnitude will require, in addi- 
tion to another order of magnitude improvement in r^t , a 
significant experimental improvement in the muon mass 
determination and a significant improvement in the var- 
ious theoretical corrections. 
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This cancellation between final-state interaction effects and bind- 
ing energy effects is familiar from work by Huff |49| and Uberall 
et al. 1501 for muonic atoms. So-called Huff factors, that repre- 



sent the difference between the ii~ decay rate in free space and 
the /i~ decay rate in muonic atoms, have been calculated across 
the periodic table. The calculated effects for bound negative 
muons are supported by experimental data from muonic atoms. 
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